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Monte Carlo Matrix Calculation with 
Punched Card Machines 


ForsyTHE & LIEBLER have described matrix inversion by a Monte Carlo 
method [MTAC, v. 4, p. 127]. While they state that the method is “best 
suited for a human computer with a table of random digits and no calcu- 
lating machine,” the principle has been adapted for use with punched card 
machines. The method has been extended to finding powers of a matrix. 
While an IBM 405 tabulator was used in this work, there is no reason why 
a newer machine cannot be used. 

The class of matrices B = (b;;) that may be inverted by this method 
is limited to those in which 
(1) the largest latent root of A* = (a;;*) is less than unity, where a;;* is the 

absolute value of a;; = 6;; — bi; and 
(2) the sum of the elements in any row of A* is less than unity 
In certain cases these restrictions may be removed. 

In this modification of the method of Forsythe and Liebler, the illus- 
trative urns containing numbered balls are replaced by a group of punched 
| cards. Each card is divided into two principal fields: a randomizing field 

consisting of 30 to 50 random digits and a playing field prepared for each 

matrix as described below. Each of the card columns of the playing field 

of the deck corresponds to one of the urns. Each of the m different digits 

punched in the card column corresponds to the numbered balls in that urn. 

An x-punch in a card column corresponds to a “stop ball” in the urns of 

‘the original description. The deck of cards is treated as a unit and is so 
prepared that the probability of a 7 appearing in the tth card column is p,;. 

In the case where A has both positive and negative elements, a value of —1 

is attached to each draw corresponding to a negative element, +1 to all 

others. The negative value is indicated by a 12 punch appearing in the 

corresponding card column above the row punch. 

In terms of the original reference, m games are played simultaneously 

to yield one row of the desired inverse B-'. The passage of each card through 

| the machine constitutes one draw in the game.' The machine is wired so 
that the original gambling procedure is followed. Instructed by a digit 

emitter, the machine reads the card column corresponding to the ith row. 

The digit punched in that column activates the mechanism which selects 

one column of the next card fer reading. The digit in each column selected 

causes the reading of the corresponding column in the following card. This 

| procedure is continued until an x is read. The column from which the x 
was read scores +1, all others score zero. As the values are to be multiplied 

together, the score is —1 if an odd number of negative draws was made, 

otherwise the score is +1. A designated counter keeps count of the total 

games played. The machine does not stop after each game, but starts anew 

on the card following the x card, again instructed by the digit emitter. 

After the entire deck of cards passes through, a total is taken and the 

machine prints total draws, total games and net score for each of the a 
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games. The actual elements of the inverse matrix are obtained by dividing 
by the product of the stop probability times the number of games played. 
If the stop probability can be taken as 0.1, and the number of games as an 
exact power of 10, the machine will actually list the elements found. By 
playing successively, starting with ¢ = 1 and continuing through i = n, the 
entire inverse will be listed with the elements in their correct positions. 

By a simple modification, powers of certain matrices may be found by 
playing a similar Monte Carlo game. In this case, the stop probabilities 
are no longer used to terminate each game. Instead, exactly r + 1 draws 
are made and the scoring is performed as before. As all possible routes 
from 4 to j consisting of exactly r steps will be covered in a large number of 
games, this is tantamount to raising the matrix to the rth power. The 
restriction on the size of the largest latent root no longer applies here, but 

n 


it is required that }>a;; = 1. The instruction to score every (r + 1)st draw 


j=l 
may be given the machine by interspersing control cards by means of a 
counting collator or by sequencing using interlocked selector relays as 
described below. 

The order of matrices that may be raised to powers or inverted depends 
on the selector capacity of the machine. For general functions, one single- 
position selector, one three-position selector and one ten-position selector 
are required as well as two digit-selectors. For matrices of the mth order, 
n two-position selectors and m single-position selectors are also needed. 

The advantages inherent in this method are: 


1. It may be used with a common machine present in nearly all punched 
card equipped installations. 

2. The actual working procedure is far simpler than for any previously 
reported matrix inversion method. 

3. Approximate inverses may be obtained rapidly. The accuracy of a 
good approximation to an inverse may be easily improved.? 

4. Using a machine with large selector capacity,* approximate inverses 
and powers may be obtained for high order matrices. In going from a 
matrix of the mth order to one of the (2”)th order, twice as many cards 
should go through the machine in twice the number of passes. Thus, as 
originally stated, the number of operations is proportional to the second 
instead of the third power of the order, as is otherwise the case. 

5. The process of finding all powers from 2 to r is an exceedingly simple 
one. It should be useful in evaluation of latent roots and vectors. 


Example 1—Inversion of a seventh order symmetrical matrix 





of —.1 1 —.1 —.1 —.1 —.1 

-.J 7 —.1 —.1 —.1 -.1 —.1 

-.1 -.1 7 J —.1 -.1 —.1 
B={|-.1 —.1 1 7 -.1 —.1 -.1 
—.1 J —.1 —§ 7 —.1 —.1 

—.1 ~.1 =~ .2 =-.1 -.1 7 -.1 

-.!1 —.1 -.1 -.1 -.1 —.1 P 
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ing Be Al mt = a 1 | A 
ed. A=| A 3 A A A A Al pj=] 1 
an 1 1 3 1 a 1 | a 
By etc etc. 
the 
Bo" =(|2.67 1.33 1.66 1.33 1.55 1.00 1.33] 

by etc. 
ties (90 games per element) 
1W 
a Bo =|2.536 1.224 1.288 1.275 1.263 1.263 1.263] 
- of etc. 
“he (784 games per element) 
m Bo =|2.500 1.250 1.250 1.250 1.250 1.250 1.250| 
‘aw etc. 
' (calculated values) 
fa 
- Example 2—A seventh order symmetrical matrix raised to the seventh power. 
ids 88 02 02 02 02 02 .02 
sle- 02 88 02 02 02 02 .02 
tor .02 02 88 .02 .02 .02 .02 
ler, A =|.02 02 02 .88 02 .02 .02 

02 02 02 02 88 02 02 

02 02 02 .02 02 88 02 

4 02 02 02 .02 02 02 88 

e 

440 .120 .088 072 .080 .088 112 
sly A’ =|.121 A27 096  ~—-«.081 081 .096 .096 

11 .079 428 -103 .095 .087 .095 





fa etc. 
(125 games per element) 
ses 


ta A’ =| .4388 .0933 .0954 .0952 .0946 .0946 .0946 | 
rds etc. 

as (825 games per element) 

d 
. A’ =| .4410 .0932 .0932 .0932 .0932 .0932 .0932| 
ple etc. 


(calculated values) 


Example 3—A seventh order nonsymmetrical matrix raised to the sixth 
power. 


88 02 .02 02 .02 —.02 .02 
02 88 02 002 02 —.02 .02 
02 02 88 .02 02 —.02 .02 
A=; 02 .02 .02 .88 .02 —.02 .02 

02 02 02 .02 88 —.02 .02 

02 02 02 02 .02 —.02  .88 
—.02 —.02 —.02 —.02 —.02 .88 —.02 
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At= | 4683 .0739 .0880 .0704 .0704 —.0422 —.0523| 
(568 games per element) 


A782 .0736 §=©.0736)3=— 0736 = .0736 —.0025 —.0025 
0736 4782 .0736 .0736 .0736 —.0025 —.0025 


AS= 








~.0025 —.0025 —.0025 —.0025 —.0025 —.0025 .4782 
(calculated values) 


Details of the procedure.—Card Preparation. 

A. Random numbers may be keypunched from standard tables‘ or 
pseudo-random numbers may be produced [MTAC, v. 5, p. 4]. 

B. Preparation of a group of cards correctly punched in the playing 
field may be accomplished as follows. 


1. Prepare the needed table of probabilities by listing the absolute values 
of the matrix with rows and columns transposed. (To prepare cards for 
matrix inversion, include the stop probabilities as the final row in the table.) 
The probabilities in each column (including stop probabilities if present) 
should add to unity. 

2. An effective, but tedious, procedure for preparing the playing field 
for a single card column would be to partition the deck by actual count 
into m (or m + 1 in the case of inversion) groups each containing ,;% of 
the card count. For each group the appropriate value of 7 would be gang 
punched into the card column representing i. A 12-punch is included when 
a value of —1 is attached to the probability and an x-punch represents 
the stop cards. 

3. This may also be done by preparing a table of accumulated proba- 
bilities, changing this to a table of turning points and collating with cards 
sorted on a d digit random field when the probabilities are given to d signifi- 
cant figures. The values of j as well as any x or 12 punches are interspersed 
gang-punched onto the playing field of the random number cards. 


Example of Preparation of Table of Turning Points for One Column in the 
Inversion of a Fourth Order Matrix 





Column of Accumulated Turning 
j Probabilities Probabilities Points j 
1 .323 .323 000 1 
2 — .076 .399 .323 2* 
3 458 .857 .399 3 
4 .109 .966 .857 4 
Pi .034 1.000 .966 x 
Keypunch 


* (also 12 punch) 
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4. Using either of these procedures, the cards should be sorted on several 


random fields between the preparation of successive columns of the playing 
field. 


Performing the Calculation 


Place the designated wire in position 1 of digit selector 2. Pass entire 
deck of cards through the tabulator. Press hand key for total. Successively 
pass cards through with wire in digit selector positions 2 to n. 


Machine Wiring 


1. Counters. Use n + 2 counter groups. Read card count into all groups. 
Add and subtract hubs go to selector F (NX and X, respectively), from C 
of that selector, through the delayed selectors and the entry selector to 
“plug to C.’’ Net balance and clear all counters on final total, emitting an 
asterisk next to negative totals. One counter group is for card count; another 
counts games played (plug to C through entry selector). 

2. Reading Selectors. n two-position selectors are used. Each X goes to 
its corresponding upper brush position in the playing field. C is common 
and leads to the entry selector. The X in the other position goes to the 
D-PU of the corresponding delayed selector. The C is again common and 
leads to a digit of DS 2. (UCI goes to DSC 2, making DS 2 a digit emitter.) 
The PU of each selector is through the corresponding hub of DS 1. 

3. Delayed Selectors. These are picked up by a digit impulse passing 
into the C of all reading selectors and through the X positions to the D-PU 
of the corresponding selector. Thus, delayed selector 3 picks up one cycle 
after reading selector 3 picks up. (This is done so that the machine can 
store the number of the “winning” digit.) The delayed selectors carry a 
common ‘‘plug to C” on their C hubs: the X hubs lead to the corresponding 
C positions in selector F. 

4. Entry Selector. Three-positions picked up by the first card or by an 
11 from Digit Selector 1. One position carries a ‘“‘plug to C’’ impulse to the 
counter which records the number of games played and to the C of the 
delayed selectors. A second is used as a drop-out for the negative impulse 
counter. The third selects either the output of the reading selector bank-NX 
or an emitted digit from DS 2 (which is impulsed by UCI)-X to DSC 1. 

5. Negative Impulse Counter. To distinguish an odd from an even number 
of 12 impulses, the selector picks up on the first negative impulse, drops 
out on the second, picks up on the third, etc. An 11 emitted by DS 2 passes 
through the entry selector (C to NX), through a special drop-out selector 
(C to NX), through one position of Selector F (C to X) to the X—PU of 
Selector F. The 12 position in DS 1 is wired through Selector F (C to N) 
to the D—PU of that selector. The corresponding X position carries the 12 
impulse to the special drop-out selector. By this means, Selector F remains 
picked up after odd 12’s, and drops out on all even 12’s and as each game 
starts. 

6. Sequencing Circuit. By passing an x through one distributor (C to X) 
to the X—PU of the next distributor, a chain of r + 1 distributors may be 
linked together. If the (r + 1)st distributor picks up the entry selector and 
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starts the chain over again, every (r + 1)st draw will be scored. This circuit 
is used for raising matrices to powers without employing control cards. 


ASCHER OPLER 
Dow Chemical Co. 


Pittsburg, Calif. 

1 This is actually sampling without replacement and therefore does not strictly conform 
to the MARKov process originally described. However, if the ratio of the number of punched 
cards to the order of the matrix is very large, nonreplacement has negligible effect on the 
remaining transition probabilities until the exhaustion of the deck approaches. 


H. Hote tino, “Some new methods in matrix calculation,”” Ann. Math. Stat. v. 14, 
1943, p. 1-34. 
3 Phe IBM 101 statistical machine should be able to invert matrices up to the 30th 
order because of its selector capacity. 


a Random numbers on punched cards are available from the RAND Corp., Santa Monica, 
lif. 


Computing Logical Truth with the California 
Digital Computer 


1. Introduction. A problem which occasionally arises and in quite sur- 
prisingly diverse fields of endeavor is the computation of a truth table for 
a sentence built up of simple sentences connected by the simple sentential 
connectives. The problem is to find the truth value of the complete sentence 
for all combinations of truth values for the component sentences. If the 
complete sentence is a combination of ” two-valued sentences, then it itself 
is two-valued and 2* possible conditions must be considered. Thus if the 
number of different component sentences -is at all large, an extremely 
tedious computation is necessary. 

Since the computation itself is performed in a routine manner, the 
possibility of using a high-speed, automatic computer to carry through the 
details suggests itself. Indeed, a machine for this purpose (the Kalin- 
Burkhart Logical-Truth Calculator) has already been built. It seems 
desirable, however, to investigate the possibility of solving these problems 
using general-purpose, digital computers as many of these will be in opera- 
tion within the next few years. With this in mind, a program for truth-table 
calculation is worked out in the following for the California Digital Com- 
puter [MTAC, v. 5, p. 57-61]. The general plan should be applicable to 
any digital computer. 

Most problem-solving programs require three basic routines. First the 
solution of the problem with a given set of data, second storing the result 
at the proper location in the memory, and third altering the data as required 
for the next computation. These will be discussed in subsequent sections. 

Since the details of a computation program depend to a considerable 
extent on the characteristics of the computer for which it is developed, a 
brief description of the operations needed and available in the California 
Digital Computer (commonly called Caldic) will be presented. Other opera- 
tions than those described are also available, but the discussion will not be 
complicated by describing them. 

The program and data are given to the machine in the form of holes 
punched in a standard teletype tape. Results are printed out of the machine 
on a similar tape. Numerical data and results are in the form of decimal 
numbers which are used by the machine in a binary-coded form. The fact 
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that numbers available to the operator are decimal is of importance in 
programming problems. This is particularly true of problems of the sort 
considered here which are essentially binary in nature. The numbers contain 
ten decimal digits and algebraic sign. 

Before computation starts, data and orders are all transferred from the 
tape to the internal memory of the machine. This transfer takes approxi- 
mately 18 milliseconds per number. The internal memory is a magnetic 
drum with a capacity of 10000, ten digit numbers. Each memory position 
is identified by a four digit, decimal number. No differentiation is made 
between orders and numbers as far as the memory is concerned. The number 
stored in a memory box is not altered by any operation except that of 
putting another number in the box. 

The first two digits in an order identify the operation to be performed. 
(The operator sees and uses them as letters, the machine as numbers.) The 
next four digits usually are the address of the operand in the memory. 
Exceptions will be discussed with the particular operations concerned. The 
final four digits are zeros. 

The results of computations appear in the A register. It takes a special 
order to transfer them to the memory if this is necessary before another 
operation is to be performed. 

The other register of interest is the order counter. The last entry on the 
input tape is a four digit number which is stored in the order counter and 
is the address of the first order. Thereafter as each order is completed the 
number in the order counter is increased by one to give the address of 
the next order. An exception is when a special order is used to change the 
number in the order counter. 

The operations of interest are described in the following. The time for 
the operation should be understood to be the average time for a large number 
of random repetitions of the operation. 


1. ad M. Average time 17.6 milliseconds. M is a four digit number. 
Add the number in the memory box M to the number in the A register. 
This operation takes proper account of the algebraic signs of the two 
numbers to be added. In case the operation causes the A register to over- 
flow, the machine stops unless the next order is the cp order. (See below.) 
At the completion of the ad operation the sum stands in the A register but 
without the overflowing digit if such exists. 

2. su M. Average time 17.6 milliseconds. Subtract the number in the 
memory box M from the number in the A register. At the end of the opera- 
tion the difference stands in the A register. The remarks about sign and 
overflow under ad apply here also. 

3. mh M. Average time 32.4 milliseconds. Multiply the number in 
memory box M by the number in the A register. This operation takes proper 
account of the algebraic sign of the two factors. It proceeds as if the decimal 
point is just to the left of each of the two factors. Therefore no overflow 
can occur. 

4. tm M. Average time 17.3 milliseconds. Transfer the number in the 
A register to the memory box M. This operation leaves the A register clear. 

5. sln. Average time 9.1 milliseconds. n is a four digit number with 
the first three digits always zero. Shift the number in the A register n 
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columns to the left. The first n digits after the sign are lost and n zeros are 
added at the right. 

6. srn. Average time 8.8 milliseconds. n same as in the sl operation. 
Shift the number in the A register n columns to the right. The n right hand 
digits are lost and n zeros are inserted after the sign. 

7. cpn. Average time 8.8 milliseconds. n is any four digit number. If 
an overflow occurred during the preceding operation, replace the number 
in the order counter by n. Otherwise ignore this order. 

8. spn. Average time 8.8 milliseconds. n is any four digit number. 
Replace the number in the order counter by n. 

9. ca M. Average time 17.5 milliseconds. Clear the address part of the 
order in the memory box M and add the remainder of the order to the 
number in the A register. 

10. po M. Average time 1117.2 milliseconds. Print out the number in 
memory box M. 


In general some of the above orders would also involve an additional 
register, the R register. However, in all the computations discussed in this 
paper the R register will hold nothing but zeros so may be ignored. 

2. Problem Solution. In the following discussion marks used to dis- 
tinguish a thing from its name will be omitted unless necessary to avoid 
ambiguities. 

The operations on or between sentences to be considered are: negation, 
conjunction, disjunction, material implication, and equivalence. For pur- 
poses of the following discussion the word, combination, will represent any 
one of the latter four. 

It will be convenient to make use of the following conventional terms in 
describing the problem solution. A schema is the result of writing in a row 
the following symbols in any order: (,), »;, and the signs for negation and 
the various combinations. The sign, p;, represents the various component 
sentences, j taking on various integral values to distinguish between them. 
Only well-formed schemas, i.e. schemas in which the parentheses are used 
to give meaningful groupings, will be used. 

The various ~;’s are, by definition, schemas of zero order. A schema 
of order n + 1 is constructed either by putting a negation sign before a 
schema of order n (which itself must be set off by parentheses) or by forming 
a combination of two schemas (each set off by parentheses) one of which 
is of order n and the other of which is of order not greater than n. 

It is assumed that the truth values of all schemas of zero order are 
available. Computation proceeds first by computing the truth values of 
schemas of order one. Then from the values of schemas of order zero and 
of order one the truth values of all schemas of order two.can be calculated. 
This process continues until the complete schema has been evaluated. The 
result is transferred to the memory, the ,’s are altered to the next set of 
values to be used, and the process is repeated. 

This may be illustrated by the following example. Suppose the fifth 
order schema 


((p1 — Po) A (((do A pi) A ~ Ps) — P2)) A ((p2 > pi) A (~ bo > ~s)) 


is to be evaluated. (Parentheses around individual variables and their 
negations are omitted.) Let a; represent schemas of order one, b; schemas 
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of order two, and so on. Then 


a, + (pi: Po) b, <> (ds A ae) 
a, +> (~)s) bz <> (a4—-a2) 
az +> (p2—> f) C1 > (bz A as) 
a4 +> (~ po) C2 +> (b;— pr) 
a,< (Po A Pi) d, o> (Ce A a) 


e+ (d, AG) 


Within any order-group any order of calculation may be used. The order 
shown is convenient since ds, bz, ¢, and d; are used only in the succeeding 
computation and thus need not be transferred to the memory. This possi- 
bility should always be exploited since in complicated problems it can lead 
to appreciable savings in computation time. 

It remains to work out the evaluation of the negation and combinations. 
Many procedures have been tried, and those tabulated below require less 
computation time than any others considered. The times given include 
transfer from and to memory which may sometimes be avoided as discussed 
above. The operations to be performed follow from left to right in sequence. 
The symbols p and q are used for the component sentences. The results 
of the operations for the various cases are tabulated below the operation. 
The symbol (x) indicates ‘‘the address in the memory of the number x.” 
The symbol C(M) indicates ‘“‘the contents of the memory box M.” The 
symbol 0*1 indicates ‘‘x zero’s followed by the digit 1.” 


Negation (52.5 milliseconds) 
Pp ad(01) su(p) tmM C(M) 
00 = «(Ol 01 00 01 
01 O01 00 00 00 


Conjunction (76.4 milliseconds) 


p q ad(p) mh(q) sl0002 tmM C(M) 

00 «600 00 0000 00 00 00 

00 «(Ol 00 0000 00 00 00 

01 00 01 0000 00 00 00 

01 ol 01 0001 O01 00 01 

Disjunction (88.0 milliseconds) 

p q ad(p) ad(q) ad(09) sr0009 sl0008 tmM C(M) 

00 «600 00 00 09 0°0 00 00 00 

00 = «(Ol 00 01 10 0°1 01 00 01 

01 00 01 01 10 0°1 01 00 01 

01 Ol 01 02 il 0°1 01 00 01 

Material Implication (88.0 milliseconds) 

p q ad(q) su(p) ad(10) sr0009 si0008 tmM C(M) 

00 «600 00 00 10 0*1 01 00 01 

00 «Ol 01 01 11 0°1 01 00 01 

01 00 00 —01 09 0°1 00 00 00 

01 Ol 01 00 10 0°1 01 00 01 

Equivalence (109.8 milliseconds) 

Pp q ad(p) ad(q) ad(98) cp0008 ad(01) cp0009 sl0002 ad(01) tmM C(M) 
00 «62400 00 00 98 98 99 99 00 01 00 01 
00 = «=O 00 01 99 99 1.00 00 00 00 
01 00 01 01 99 99 1.00 00 00 00 
01 Ol 01 02 1.00 00 01 


00 01 
Order counter 0001 0002 0003 0004 0005 0006 0007 0008 0009 
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It will be noticed that 01 is the truth value of a true sentence and 00 is the 
truth value of a false sentence. The zero in the first place simplifies the com- 
putation of the implication and has no effect on other computations. 

3. Handling of Results. Since results consist of a single digit, it is possible 
to store ten of them in one memory box. Thus the first answer may be put 
in the left hand digit position, the second in the next to the left and so on 
until the memory box is filled. Then the contents of the box may be printed 
out, the box cleared, and the process repeated. 

Storing ten answers in each box reduces the number of print operations 
by a factor-of ten which is desirable since printing is by far the slowest of 
the machine operations. 

Printing during computations makes the results available during com- 
putation which may be convenient. Also if the machine breaks down during 
computation, it will not be necessary to recalculate completely. 

4. Altering Data. In order to make use of the advantages of automatic 
machine calculation, the machine itself should compute all but the initial 
data. In this problem the data are in the form of (binary) truth values and 
it is desirable to carry out the computation for all possible combinations 
of truth values. 

Let the truth value (0 or 1) of each component sentence be p; where 
j ranges over the values from 0 to m — 1, being the number of component 
sentences to be considered. Construct a binary number k the left hand 
(or most significant) digit being the result of writing the value of p,_1, the 
next most significant digit being the result of writing the value of p,_2, and 
so on. The least significant digit will be the result of writing the value of po. 
All combinations are covered if & starts as zero and then is increased by 
unity for each computation until every digit of k is one. 

In machines using binary arithmetic this procedure can be worked out 
quite readily. The machine being considered here uses decimal arithmetic 
which complicates the problem somewhat. The procedure finally adopted 
makes use of the theorem? which follows. For economy in writing, ‘‘p,’’ will 
be considered equivalent to ‘‘the result of writing the value of p;.”’ 

Definition. A binary number b is defined as the array 


Pn-1Pn—-2 *** Pipo. 


Theorem. The number is replaced by k + 1 if and only if for every 7 
(if0 <7 < m —1, then 9; is replaced by 1 — 9; if and only if for every m 
(if 0 < m <j then p,, = 1)). Or more compactly 


R(B)k + 149;,(0 Sj Sm — 1->(p,(B)1 — pjon(0 S m <j Pm = 1))) 


(It is understood that p; remains unchanged unless it is replaced by 1 — 9; 
by the above rule, and that at every step fp» is replaced by 1 — pp.) 

This procedure is programmed quite readily. 

5. Terminating the Problem. The machine will not stop, of course, 
unless it has definite instructions to do so. These may be obtained readily 
by adding another digit », to the left of k. This is zero throughout the 
computation but will change to one when a one is added to the last value 
of k to be considered. Thus in the preceding section when the largest value 
of j satisfying the last equivalence equals m, the computation is finished and 
the machine may be instructed to stop. 
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1e The situation may occur in which all possible combinations of truth 
n- values need not be considered. For instance, only true or only false values 

of a particular sentence may be of interest. In this case this value will be 
le held as a constant (and not stored in the memory with the other p’s) and 
ut n will be one less than the number of component sentences. 
on 6. The Program. The operations discussed in the four preceding sections 
ed are embodied in the flow chart, figure 1, and together make up the program. 
When the flow chart is worked out in sufficient detail, the program is readily 
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coded or translated into a language to which the machine will respond. 
A code for this problem is given in Appendix I. 

Within the rectangles of the flow chart are first the memory boxes in 
which the instructions to perform the operation demanded in the rectangle 
are to be stored (for the particular code given in Appendix I). Secondly 
the operation to be performed is given. In case the operation to be per- 
formed is in the form of a question, the alternative answers to the question 
indicate the paths by which the computation will proceed. The capital 
letter in the rectangle identifies the operation for the discussion below. 

The expression x(B)y means that at this point y is to be substituted for x 
wherever x appears in the program. It may be read ‘‘x becomes y.”’ 

Operations D and E must be performed a number of times equal to the 
largest digit which is less than 0.1 X 2". Operations H and I must be per- 
formed a number of times equal to 2"*! — (m + 2). Operation L must be 
performed once. All other operations must be performed 2* times. 

Operation A is self-explanatory. Operations B, C, D, and E are for 
handling the results. Operations F, G, H, I, and J take care of altering the 
data. Operations K and L are for terminating the problem. 

7. Conclusion. The program given here was developed and coded with 
the idea in mind of minimizing the changes necessary in going from one 
problem to another. It was thought that 19 variables is beyond the limits 
of practicality (for this machine, at least, since 19 variables would take 
over six days of machine time). As long as no more than 19 variables are 
used the instructions in memory boxes 0000 to 0072 inclusive will apply 
to any problem. Memory box 0073 will contain the number 1.00—0.0in 
(decimal point to the left of the most significant figure of numbers in the 
memory). This will be used to terminate the computation. Starting with 
memory box 0074 the schema is evaluated as discussed in section 2. The 
schema given there as an example is the one coded in Appendix I. The 
truth value of the schema is finally left in the left hand digit position of 
the A register. (Note that the routines of section 2 leave the value in the 
position next to the left hand one so that the last operation must be altered 
accordingly.) Then the command sp0021 is used to go to operation B of 
figure 1. 

In preparation for working these problems a master tape going only 
through 0072 will be prepared. For a particular problem this will be dupli- 
cated and the additional orders needed will be added on the duplicate. 
Then the entire contents of the tape will be fed into the machine, the A and 
R registers cleared, and the machine started with the order counter at 0074. 

Using the average times given in section I the problem time is 


T = 0.91 X 2” — 0.18” + 0.75 + 0.077c X 2" seconds, 


where c is the number of connectives in the schema. It should be empha- 
sized that this is a rough estimate and depends on the assumption that on 
the average it will take half a drum revolution to find the memory box 
requested. In the example given here m = 4 and c = 11 (~ 3 need be 
calculated only once) yielding T = 26 seconds. (It is understood, of course, 
that this short example would be worked more economically manually.) 
Available data on the Kalin-Burkhart machine are rather sparse, but 
it would be interesting to compare its operation with that of the program 
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developed here. It would appear that the Kalin-Burkhart machine will 
compute the problem a little faster than will the Caldic. The Caldic, how- 
ever, is rather slow as electronic computers go. Problem preparation time 
should be about the same. On the other hand the Kalin-Burkhart machine 
is limited to problems with 12 component sentences and 11 connectives. 
With the Caldic the number of component sentences is limited only by the 
time for which the machine is available for the particular problem and there 
appears to be no practical limitation on the number of connectives. 

Appendix I. The Coded Program. These quantities will be fed in sequence 
into the memory the first into box 0000, the second into box 0001, and so on 
until the list is exhausted. It is understood that the machine will treat 
letters as digits and that each number will be filled out with zeros to make 
ten digits. Every quantity has a plus sign attached. Under each memory 
box number is given its contents. 


0000 0001 0002 0003 0004 0005 0006 0007 

0 0 0 0 0 0 0 0 

0008 0009 0010 0011 0012 0013 0014 0015 

0 0 0 0 0 0 0 0 

0016 0017 0018 0019 0020 0021 0022 0023 

0 0 0 0 0 sr0000 ad0020 tm0020 
0024 0025 0026 0027 0028 0029 0030 0031 


ad0021 sl0005 ad0069 cp0032 sr0005 ca0021 tm0021 sp0036 


0032 0033 0034 0035 0036 0037 0038 0039 
po0020 tm0020 ca0021 tm0021 ca0043 tm0042 ca0047 tm0047 


0040 0041 0042 0043 0044 0045 0046 0047 
ca0056 tm0056 ad0000 ad0071 cp0056 s10002 ad0067 tm0000 


0048 0049 0050 0051 0052 0053 0054 0055 
ad0042 sl0004 ad0073 cp0054 310002 sp0074 po0020 STOP 
0056 0057 0058 0059 0060 0061 0062 0063 
tm0000 ad0042 ad0072 tm0042 ad0047 ad0072 tm0047 ad0056 
0064 0065 0066 0067 0068 0069 0070 0071 
ad0072 tm0056 sp0042 01 09 10 98 99 

0072 0073 0074 0075 0076 0077 0078 0079 
000001 96 ad0000 su0001 ad0069 sr0009 310008 tm0116 
0080 0081 0082 0083 0084 0085 0086 0087 
ad0067 su0003 tm0117 ad0001 su0002 ad0069 sr0009 sl0008 
0088 0089 0090 0091 0092 0093 0094 0095 
tm0118 ad0067 su0000 tm0119 ad0000 mh0001 sl0002 mh0117 
0096 0097 0098 0099 0100 0101 0102 0103 
sl0002 tm0120 ad0117 su0119 ad0069 sr0009 sl10008 mh0118 
0104 0105 0106 0107 0108 0109 0110 0111 


10002 tm0121 ad0002 su0120 ad0069 sr0009 310008 mh0116 


0112 0113 0114 0115 
10002 mh0121 s10003 sp0021 
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Appendix II. Glossary. 


nm the number of component sentences in the complete schema. 
p; the jth component sentence. 

a;, b;, etc. the jth component schema of order 1, 2, etc. 

c the number of connectives in the schema. 

T estimated computing time in seconds. 

1 (or 01) in the proper context the truth value of a true sentence. 
0 (or 00) in the proper context the truth value of a false sentence. 
x(B)y y is to replace x wherever x appears. 

(h) address in the memory of the number h. 

C(M) contents of the memory box M. 

~ Itis not the case that... 


A _ ea 
ee 
— If... then. 
o . if and only if... 


« Foreveryx... 
WiLTon R. ABBOTT 
Univ. of California 
Berkeley 


1 Epmunp C. BERKELEY, Giant Brains, New York, 1949, Ch. 9. 

2 GEORGE W. PatTERsON, Logical Syntax and Transformation Rules. Moore School of 
Electrical Engineering, Research Division Report 50-8, Univ. of Pennsylvania, Phila- 
delphia, 1949. 


On the Accuracy of Runge-Kutta’s Method 


1. Introduction. While the accuracy of the most frequently used methods 
of integrating differential equations is fairly well known, that of the Runge- 
Kutta method does not seem to be too well established; except for a formula 
in Bieberbach’s text! on differential equations there are no references per- 
taining to the error inherent in the Runge-Kutta method to be found in the 
standard textbooks on this subject. 

Since this method may be employed quite advantageously in many cases 
of practical interest it is important to have on hand an estimate of the error. 
The purpose of the following sections is to provide such an estimate. As a 
comparison shows, the bound derived for this error seems to be somewhat 
better than the one cited by Bieberbach. 

2. Runge-Kutta’s Fourth Order Method. In trying to find that solution 
of the differential equation 


(1) dy/dx = f(x,y), — y(xo) = ¥, 


at x, = x9 + h, which agrees with the exact Taylor expansion about xo: 


(2) y(x1) = yo + hyo’ + h*(y0"/2) + hP(y0'”/6) + h*(yo'"/24) 
+ h5(yo"/120) + ++: 
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up to the term in hk‘, Runge and Kutta developed the following formulae: 


(3) y(x1) ~ Yo + (hi + 2ke + 2h + ha) /6, 

ki=hfo, fo = f(x,y), 

ky = hf (xo + h/2, Yo + k,/2), 

ks = hf(xo + h/2, vo + k/2), 

ke = hf(xo + h, yo + ks). 
To get an estimate of the truncation error inherent in this procedure, one 
may apply the method first to an interval of length h, = 4, and then 
integrate over two consecutive intervals of length 4, = h/2. Having the 
results Y;, Y2 of these integrations it is easy to obtain an estimate of the 


error of the second integration: Since the values Y,, Y2 differ from the exact 
value y; by certain errors E;, E: 


¥Y,=y7+ E, Y, = 9 + E,, 


where 

E,= Ch‘, E, = 2ChS = E,/16, 
obviously 
(4) E, = (Y; — ¥,)/15. 


3. Calculation of the Error Term. A more accurate estimate of the error 
is obtained by a comparison of the exact coefficient (yo"/120) of h°, as it 
occurs in the Taylor expansion (2), with the approximate one resulting from 


Runge-Kutta’s algorithm (3). Suppose equations (3) have been expressed 
in the form 


y(x1) = Yo + Crh + Crh? + Cyh® + Ci + Ch + --- 
Then C; = y/i!, i = 1, 2, 3, 4, 
(5) Cs = (w"/5!) + €, 


and the following considerations are concerned with the determination of e. 
In computing the successive total derivatives of a function u = u(x, y(x)), 
y’ = f, it is helpful to make use of the operator? 


D = 4/dx + fa/dy. 
This operator has the following properties: 


D> = ¥ G)fio"/ax"-*ay, 
D(Du) = Du + n(Df)D™"u,, 
(6) Dy = u. 


Let us now apply the operator D to f(x, y(x)) and use the abbreviating 
notations 


T* = D*f, S* = D*f,. 
Then (6) may be expressed as 
7 D(T*) = T*#! + nTS™. 
7) DS = 3 + Thy 
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Applying (7) repeatedly we get 
yes 
y” -_ Df = TT 
y" = D(T) = T?+ Th, 
y*” = D(T? + Tf,) = D(T?) + D(T)f, + TS 
= T* + f,T? + 3ST + Tf. 
Similarly, 
y = T*+ f,T? + 6TS + 4ST? + f7T? + 3f,(P)? + 7f,ST + fT, 


with (P)? = (Df)?. 
The expansions of the Runge-Kutta’s expressions (3) are somewhat 
more laborious to calculate. First we find that 


kee = hf(xo + h/2, yo + ki/2) = hLf + (h/2)fe + (hi/2) fy 
(8) + 3((h?/4) fox + 2(h/2)(hi/2) fey + (hr?/4) fa) + °°° J 
= ky + (T/2)h? + (T?/8)h® + (T*/48)h* + (T*/384)h5 + O(h'). 
Here, as well as in the following, the arguments of f and f,‘,i are xo, yo. 
Next we compute 
ks = hf(xo + h/2, yo + ko/h) = hLf + (h/2)fe + (he/2)f, 
+ 3((h?/4) fee + 2(h/2)(ke/2) fey + (hs?/4) fv) + °°* J 


= ALf + (1/2)(hfe + Refy) + (1/8) (fos + Zhhofey + kof) + --- J. 


It is convenient to introduce here the operators: 
G; = ha/dax + kd/dy, ‘= ¥ 2, 3, 
Us = Gf = (hd/dx + kd/dy)"f. 
In terms of the U;" we may write 
(9) ks = WL f + (Us/2) + (Us?/8) + (Us8/48) + (Us*/384) + +++]. 
In order to express k; in powers of hk we notice from (8) that 
koe = ky? + fTh® + (1/4)(P + fT)h* + O(h'). 
Further, 
ke = ky + (3/2) f°Th* + O(h'), 
kot = ky! +- O(h*). 
Consequently, 
U, = hT 
U2 eal hf + Rofy 
U; + (1/2) f, Th? + (1/8) f,T*h® + (1/48) f,T*h* + O(h*) 


U2 = U;? + TSh + (1/4)(ST? + Pfy,)h* + O(h') 
US = Us + (3/2) TSh* + O(h') 
U4 = Ui! oe O(h*). 


Therefore, by (9), 
(10) ks = ky + 4TH? + Resh® + Resh* + kash’ + O(h), 
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with 
kes = (27f, + T*)/8 
ku = (3f,T*? + 6ST + T*)/48 
ko = (4f,T* + 12(ST? + Phy + TS*) + T*)/384. 


Finally we compute 


ky = hf(xo + h, Yo + ks) 
(11) = hLf + hfs + kafy + 2 fee + Zhkafey + ki? fy) + -+- J 
= hLf + Us + (Us#/2) + (U#/6) + (Us*/24) + ---]. 
Since, by (10), 
ki? = (ky + (T/2)h? + Rash® + ---)? 
= ki + fTh® + (1/4)(P + fhss)h* + O(h'), 
ki = k,’, kt = ke", 
there result the following expressions: 
U; = hf: + Rsfy 
Uz + (1/4)Tf,7h? + (1/16) f,(f,T? + 2ST) he 
U? = U? + (1/2)STf,h 
U# = U, Us = Us. 


Thus it follows from (11) that 
(12) ke = ky + Th? + Rash? + Rash* + Rash® + Ok), 
with 


ll 


kas = (1/2)(Tf, + T?) 
kus = (1/6) (6hsafy + 3ST + T*) 
kus = (1/24)(24(Risfy + Shas) + 3Pfy + 6TS + T*). 


In the Runge-Kutta expression (3) the coefficient C; of h is thus found 
to be 


Cs = (Ris + 2hos + 2hss + has) /6 
= (36P fy, + 60T7S? + 72f,ST + 36ST? 
+ 12f,2T? + 8f,T? + 107*)/1152. 
It follows that 


e = Cs — (yo"/120) 
= (— 12Tf,? + 9Pfy + 3TS + 6STf, — 3ST? 


+ 3f7T? — 2f,T* + (T*/2))/1440, 
or, more explicitly, 


e = [— 12Tf/} + 3(f,? — S)fee + 3(2Tf, + 2ff? — 2fS — Phua)fev 
(13) + 3(3P + 2T ff, + PH? re Pfu) fw oo 2fyfeze + 3(fz ai Sfu)fezy 
+ 6ffefew + (T + 2f2) ff + (T*/2) 1/1440. 


4. Estimate of the Error. Let it be assumed, now, that in a certain 
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region B(x, y) containing (xo, yo) 


\ f(x, »)| < M, 
lfeyi] < L*1/MP, 


where M, L are positive constants independent of x, y. In that case clearly 


(14) 


IT*| < zy (3) Mi(L"/M*!) = M(2L)" 


IS} <5 (QMUL/M) = LL)", |P| < (2MLY. 


By (13), then 
(15) |E| = |eh5| < (73/720) ML*h'. 
To calculate a suitable value of Z one may first compute bounds 
L;x for | f.‘4|, then define quantities 
Lig = max ((Li+j,0/M)Y, (Lig sas)¥?, (MLi4j-2,2) 9, +++, 
(MHL y ,,)MH0}, 
for 1 + 7 < 4. Then one may put 
(16) L = max (Li, In, Ls, Ty). 
5. Example. Let us consider the differential equation 
y=xty, 90) =9, 


and restrict x,y to 0O<x< 0.2, 0< y < 0.2. The exact solution is 
y = e* — x — 1, so that y(0.2) = 0.02140276. For 4, = 0.2 Runge-Kutta’s 
method gives Y; = 0.0214, for 42 = 0.1 it gives Y2 = 0.02140257. Now for 
the above example we may take M = 1, L = 1, so that for A = 0.2, by 
formula (15), 

|E| = (1/10)(2/10)§ = 0.000032 ; 


actually the error is less than 0.000003. According to the expression ex- 
hibited by BrEBERBACH! 
(17) |E| < 6MNh‘| N® — 1|/|N — 1], 
where |f| < M, |fzi,i| < N/M*, and, further, AN < 1, aM < b, where 
|x — %| <a, |y — yo| < 6. With N = 1 formula (17) leads to 
|E| < 0.0096, 
which is considerably larger than the error estimate obtained by (15). 
6. Systems of Equations. For a system of differential equations 
yi = fix, Vis Ye); yi(Xo) = Vids 1 = a 2, 
the Runge-Kutta formulae become: 
ka = hf, fio = Silxo, V0; Yoo) 
ke = hfi(xo + h/2, yo + ku/2, yo + ko /2) 
(18) kis = hfi(xo + h/2, yo + Ri2/2, veo + Re2/2) 
Ria = hfi(xo + h, Yio + Ris, Yoo + hes) 
yilx1) = vo + (1/6)(Ra + Riz + 2is + Ris). 
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Proceeding in the manner described in the foregoing sections for the 
case of a single differential equation, and making use of the abbreviations 


Rr = Df, = Dg, o& = Dy p= Dy, w= [rP 
one obtains for «, = Cis — y"(xo, Yio, Yoo) /120 the expression 


14406 = (T*/2) — 2(f,T* + f.r*) + 3(T*(f,? + fev) + PU fvfe + fobs) ] 
— 3(ST? + Rr*) + 12LS(Tf, + rf.) + R(Tgy + gs) ] 
— 6[T(Sf, + of.) + r(Rf, + ef.) ] + (TS + +R?) 
+ 0Pfy + 2T thy + fez) 
— 12(T(f,) + fegy(2fy + g2)) + thf? + fs + feks) J 
— 15f,(T*g + 7*f). 


Assumptions similar to (14) now permit one to get a bound for ¢; A*. 
If, namely, near (xo,, Yi, Y2), 


lfi| < M, | aPtetef;/AxPy,%y2"| < Letetr/ Mer! 
for 0 < p+ q+7 < 4, it is found that 
(19) |E;| < (973/720) ML*h*. 


For the differential equation y” — yv = 1, y(0) = 0, y’(0) = 1, which is 
equivalent to the system y:' = ye, yo’ = 1+ 1, 9:(0) = 0, y2(0) = 1, the 
solution is y = e* — 1. Thus (0.1) = .1051709. With 4 = 0.1, R.K.’s 
method gives y(0.1) = .1051707, whence E = 2-10~’. In the region 0 < x 
< 0.1, 0O< y < 0.11, 0 < ye < 1.11; above estimate (19) asserts that 
|E;| < 1.5-10-. 

Max LorkKIn 


Ballistic Research Laboratories 
Aberdeen Proving Ground, Md. 


1L. BIEBERBACH, Theorie der Differentialgleichungen. New York, 1944, p. 54. 
2C, Runce & H. K6nicG, Vorlesungen tiber numerisches Rechnen. Berlin. 1924. 


RECENT MATHEMATICAL TABLES 
878[F].—H. CuaTLanp & H. Davenport, “Euclid’s algorithm in real quad- 
ratic fields,” Canadian Jn. Math., v. 2, 1950, p. 289-296. 
The tables in this paper are rather special and are used to prove a result 
of INKERI! that the field k(m') is not Euclidean for 
m = 193,241,313,337,457,601. 


There are six tables corresponding to these values of m. Each gives a com- 
plete period of a periodic algorithm showing certain inequalities which 
establish the non-existence of Euciid’s algorithm in each field considered. 


1K. Inxert, “Uber den Euklidischen Algorithmus in quadratischen Zahlkérpern,” 
Acad. Sci. Fennicae, Annales, v. 41, 1947, p. 5-34. 


879[F].—A. GLopEN, “‘Factorisation de nombres NV + 1,” Euclides, v. 10, 
1950, p. 157. 
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A small table gives the complete factorization of N** + 1 for N = 2(1)13, 
21,27,32 with references to authors. Two additional entries, for N = 16 and 
N = 64, might have been given: 


16% + 1 = 274177-67280421310721 
64° + 1 = 641-6700417- 18446744069414584321. 


The first of these is the historic factorization by LANDRY of the 6-th FERMAT 
number. For the second see Amer. Math. Soc. Bull., v. 36, 1936, p. 849. 
A third factorization corresponding to N = 36 


36 + 1 = 2753- 145601 - 19854979505843329 


has been announced by HoprENoTt and Kraitcuik [ Sphinx, v. 4, 1934, 
p. 47] but no really satisfactory proof of the primality of the largest factor 
has been given. 

D. H. L. 


880(F].—HELmuT Hasse, Arithmetische Bestimmung von Grundeinheit und 
Klassenzahl in zyklischen kubischen und biquadratischen Zahlkérper. Akad. 
Wissensch., Berlin, mat. nat. Kl., Abhandlungen, 1948, n. 2, Berlin, 1950. 
vii + 95 p., 29.5 X 21 cm. 14 marks. 


This paper contains tables for the fundamental units, the cyclotomic 
units and the class numbers of all cyclic cubic and all real cyclic biquadratic 
fields with conductor = 100. 

A systematic method is developed for finding these quantities generally 
in such fields; it is applicable to the case of non-real cyclic biquadratic fields. 
The work of H. Bergstrém!' is used and generalized. 

The analytic expression for the class number depends on the index of 
the system of cyclotomic units with respect to the system of all units in the 
maximum real subfield. For real fields this is the field itself. A cyclic cubic 
field is always real and has two fundamental units, a cyclic biquadratic 
field has three or one according as it is real or imaginary. It is shown that, 
for cubic fields, there exists a fundamental unit E such that —1 and the 
conjugates of E generate all the units. For real biquadratic fields the ex- 
istence of a relative fundamental unit E is shown which has relative norm 
+1 with respect to the quadratic subfield (which is always real). The 
conjugates of E, the fundamental unit of the quadratic subfield and —1 
generate a subgroup of index 1 or 2 in the complete group of units. In both 
cases E is characterized by a minimum property of its coordinates with 
respect to the integral base, just as the fundamental unit of a real quadratic 
field corresponds to the smallest positive solution of the Pellian equation. 

The integral base is chosen in a special way by using ideas of class field 
theory. This is done in order to determine it in an invariant way. This 
can be regarded as a generalization of the fact that in the quadratic case 
the discriminant determines the base. 

O. TAUSSKY 
NBS 


1H. Bercstrém, “Die Klassenzahlformel fiir reelle quadratische Zahlkérper mit 
zusammengesetzter Diskriminante als Produkt verallgemeinerter Gaussscher Summen,” 
Jn. reine angew. Math., v. 186, 1944, p. 91-115. 
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881(F].—E. H. NEVILLE, editor, The Farey Series of Order 1025, Displaying 
Solutions of the Diophantine Equation bx — ay = 1. Designed and com- 
piled by E. H. N. (Royal Society Mathematical Tables, volume 1). 
Cambridge, University Press; published for the Royal Society, 1950. 
xxxii, 405 p. 21.5 XK 28 cm. There is also a loose card, 20.2 K 27.2 cm., 
giving the complete factorization, on one side of the numbers 1-549, 
and on the other of 500-1049. £555. Cambridge University Press, 
$18.50, New York. 


This volume was accepted for publication by the British Association 
Tables Committee in 1947. In the following year the responsibility which 
the British Association had exercised since 1871 through that Committee 
was transferred to the Royal Society which set up a Committee of its own 
to take over the work [see MTAC, v. 3, p. 333-340]. The publication of 
the volume under the CUNNINGHAM bequest for the production of new 
tables in the theory of numbers is especially appropriate. 

The so-called FAREy series, ‘F,, of order m, is the ascending series of 
irreducible fractions between 0 and 1 whose denominators do not exceed n; 
the numbers 0 and 1 are included in the forms 0/1 and 1/1. Thus, F; is 


Beek PRR GE Ba EaAeBRaA Bae 


[en ae oe ee ae Oe. ee eo eee 
The fraction } is the middle term of every Farey Series; if two terms are 
equidistant from this, their sum is unity; they have the same denominator 
and this denominator is the sum of the two numerators. To emphasize 
this property of the series, as well as to economise in presentation, F; would 
be written: 

ae 2 We Br oe OP te Ee 

47) 3 @ 7: 3 8? 

fim SS-S: 3 2s |G 
Pages 1-400 of the volume before us are taken up with the presentation of 
the 319765 terms of Fio25, written in this form. Except for the last page 
there are 400 terms on each page, of 20 lines, each line containing 20 term- 
pairs, separated as two groups of 10, and the whole arranged in 8 groups 
of 50 term-pairs. The first three and the last three printed term-pairs of 
Fio25 as given are 


1 
2 
1 


0 1 1 sii 6512 1 
1 1025 1024 and 1023 1025 2 
1 1024 1023 512 513 1 


It is clear that all terms of an F,,, m <n, are contained in the series 
of terms of F,. 


If a:/b:, d2/b2, a3/b3 are any three consecutive terms of F, 


I. (a; + a3)/(b: + 03) = ae/bo. 
II. (@;/b1) — (a2/b2) = 1/(byb2) or (a2/b2) — (a3/bs) = 1/(bebs), bi, be, ds 
being relatively prime. 
III. Any term between a,/b, and a;/b; is of the form (pa,+ ga:)/(pb.+ qhi), 
where p and g are integers, and is irreducible if p, g are relatively 
prime. 
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IV. The total number of terms in F, is ®(”) + 1, where (7) is the sum 
function of Euler’s ¢(n). 


The work under review contains three appendices (which I shall call 
A, B, and C). 

Appendix A (p. 402-403) gives the Farey series of order 50 with decimal 
equivalents to 5D. 

Appendix B (p. 404) displays all the terms of Fg, in order to give the 
reader a general idea of a Farey series which may not be deduced from the 
400 pages filled with Fo25. 

Appendix C (p. 405) exhibits the “Farey integer-series’” of order 100, 
that is, the denominators, in order, of all the terms of the series in Fi. 
The page has to be read twice (once forward and once backward), in order 
to cover all terms. 

Let us now pause to consider some historical facts connected with the 
series under consideration, insofar as they have bearing on the contents 
of the work under review. 

a. In Ladies’ Diary for 1747, p. 34, Question 281, J. May, Jr., of Am- 
sterdam, offered the problem: “It is required to find (by a general theorem) 
the number of fractions of different values, each less than unity, so that 
the greatest denominator be less than 100.”’ Since reducible fractions are 
not excluded the number here would be in excess of the number of terms 
in Fy. Solutions are presented in Ladies’ Diary for 1748, p. 23, and for 
1751, p. 30. 

B. In Journal de I’ Ecole Polyt., cahier 11, v. 4, 1802, p. 364-368, Citizen 
Haros writes on the topic ‘“Tables pour évaluer une fraction ordinaire avec 
autant de décimales qu’on voudra; et pour trouver la fraction ordinaire la 
plus simple, et qui approche sensiblement d’une fraction décimale.” 

Haros here proved the general results I and II stated above. 

yi. In 1816 Henry Goodwyn (1745-1824) printed privately and dis- 
tributed some copies of his The First Centenary of a Series of Concise and 

Useful Tables of all the complete Decimal Quotients, which can arise from 
dividing a unit, or ceny whole Number less than each Divisor by all Integers 
from 1 to 1024. This publication (‘“‘specimen” he calls it) contained xiv + 18 
pages. Page i is simply the title page quoted above; p. ii is blank; on p. iii 
is an introductory statement signed by Goodwyn and dated March 5th, 
1816. The 18 pages of tables simply contain 6D equivalents of the various 
fractions with denominators successively 1(1)100. From this series of tables 
one could get no ready idea of the order of successive terms in a so-called 
Farey series. 

2. At London in 1818 Goodwyn sent out a publication on the upper 
part of the title page of which is the italics title quoted above; then follows 
immediately : To which is now added a Tabular Series of Compiete Decimal 
Quotients for all the Proper Vulgar Fractions of which when in their lowest 
terms, neither the Numerator nor the Denominator is greater than 100: with 
the equivalent vulgar fractions prefixed. This is repeated as title page (i) of 
the separately paged addition. viii + 32 p. Here we have, p. 1-14, not only 
all the terms in the first half of Fo, and their 10D equivalent values, but 
also the means for writing down at once the remaining terms of the series 
and their decimal equivalents. 
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Hence, this table offers the means for (a) reading off all the results of 
Appendix C in the table under review; (b) displaying the Farey series of 
order 100 rather than of order 64 in Appendix B; (c) reading off the Farey 
series of order 100 with decimal equivalents rather than the series of order 
50 with diminished decimal approximations in Appendix A. 

Furthermore, on page v of this 1818 publication of Goodwyn we find 
in his survey of the series: ‘every Fraction in the Series (since they are all 
in their lowest terms), is the exponent of the ratio of the sum of the numerators 
to the sum of the denominators, of the Fractions occurring next before and 
after it.” This is just one of the results found by Haros 16 years before. 
But in this connection Goodwyn does not mention the name of anyone else. 

5. It will presently appear why we must now turn back in our chrono- 
logical treatment. In Phil. Mag., v. 47, May 1816, p. 385-386, the geologist 
John Farey (1766-1826; see Dict. Nat. Biog.) contributed a letter to the 
editor. The title for the heading to this letter is ‘On a curious property of 
vulgar fractions.’’ The letter commences: “‘On examining lately, some very 
curious and elaborate Tables of ‘Complete decimal Quotients’ calculated 
by Henry Goodwyn, Esq. of Blackheath, of which he has printed a copious 
specimen, for private circulation among curious and practical calculators, 
preparatory to the printing of the whole of those useful tables, if sufficient 
encouragement, either public or individual, should appear to warrant such 
a step: I was fortunate while so doing, to deduce from them the following 
general property; viz.,”—and then follows the Haros-Goodwyn theorem, 
with numerical illustrations from Fs and Fy. 

What deductions may we now make? Farey was writing about May 
1816 and refers to Goodwyn’s privately distributed publication of the 
previous March, and also, to tables of “Complete Decimal Quotients.” 
It is obvious that by 1816 Farey must have seen both tables of the 1818 
publication of Goodwyn, to each of which the above quoted title applies, 
in order to make public his deduction (impossible from the 1816 table alone). 
We have noted that Goodwyn formulates this same result. 

Farey’s letter was reported in Soc. Philomathique de Paris, Bull. d. Sct., 
1816, p. 112, but Farey’s name is not mentioned. This communication 
attracted the attention of Cauchy, who gave a proof, and extension, in the 
same volume of the Bulletin, p. 133-135. As we have seen these results, 
I-II, were already published by Haros 14 years earlier. 

The result III was proved by Epwarp Sane (R. Soc. Edinb., Trans., 
v. 28, 1879, p. 287). The name ‘‘Farey Series,’’ which, as we have seen, has 
not the slightest justification historically, seems to have been first used by 
Sylvester, who also tabulated IV (Phil. Mag., s. 5, v. 15, 1883, p. 251f.). 
Well knowing all the facts in question Neville would have served scholarship 
better by entitling his book: The Haros Series of Order 1025, or The Goodwyn 
Series of Order 1025. 

The table in clockmaker Achille Brocot’s Calcul des Rouages par A pproxt- 
mation. Nouvelle Méthode, Paris, 1862, p. 49-89, is a display of Fioo, with 
each term approximated to 10D; the first half is identical with Goodwyn’s 
1818 table. (72) 

In Hutte, Hilfstafeln zur I. Verwandlung von echten Briichen in Dezimal- 
briiche . . . Third newly revised edition, Berlin, 1922, p. 24-61, J. T. 
Peters checked and extended ‘“‘Brocot’s’’ table from 10D to 11D. 
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In his Manual of Gear Design, Section One, New York, 1935, p. 148-169, 
EARLE BUCKINGHAM proposed to give a table of Fin to 8D, but there are 
numerous omissions and errors; see MTAC, v. 1, 1943, p. 92. 

On the other hand, the terms of Fix in R. M. Pace, 14000 Gear Ratios, 
New York, 1942, are complete and the corresponding 11D values more 
reliable. [See MTAC, v. 1, 1943, p. 21-23; 1944, p. 326-329, 430.] 

Until the publication of Neville’s volume the most extensive F, to be 
published was that for m = 120. The statement by D.H. L. in his Guide 
to Tables in the Theory of Numbers, Washington, 1941, p. 8, that Goodwyn 
gave Fo is therefore incorrect; see MTAC, v. 1, p. 372, 1945. 

All solutions (x, y) of the equation 


(1) bx — ay =¢, (a, b, c, integers, a, b, coprimes) 
are given by the formulas 

x = ka+ Xe, y = kb+ Ye, k any integer, 
if 
(2) bX —aY = 1. 


Hence, a/b and X/Y are consecutive terms of a Farey series, and solutions 
(X, Y) of (2), and hence of (1) are known if the neighbours of a/b in F; are 
known. Neville’s Introduction explains, with numerical illustrations, not 
only the processes by which solutions of any equation (1) may be obtained, 
however much 6b exceeds n, but also how to find desired terms in Fo. 
A “Decimal Index” table on p. xxviii—xxix, and the loose card, are useful in 
this connection. 

Solutions of (1) have important applications in different parts of the 
theory of numbers. 

The two earlier tables for finding solutions of (2) are by A. L. CRELLE 
(in his Journal, v. 42, 1851, p. 304-313), and by A. J. CUNNINGHAM (in his 
Quadratic and Linear Tables, London, 1927, p. 134-157). Crelle gives imme- 
diate solutions of (2) for a < 6 < 120. The results of Cunningham’s table, 
for b < 100, a < 100, are therefore included in Crelle’s. Within its range 
the Crelle table would naturally furnish results more readily than Neville’s 
tables would. 

To ensure accuracy in this splendid new volume of the British Mathe- 
matical Tables Committee more than 40 people collaborated with Professor 
Neville. Hence we must wonder why neither they nor the printer noticed 
that on page 399 are ten numbers which should have been in ordinary 
rather than in blackface type. 

This volume is dedicated ‘‘To the memory of SRinIvASA RAMANUJAN,”” 
Neville’s friend, since ‘‘to every mathematician of our time Farey’s series” 
recalls his name. 


R. C. ARCHIBALD 
Brown University 
Providence, R. I. 


1 Not only the publications of Goopwyn to which we refer, but others, have been dis- 
cussed in thorough fashion by ‘ W. L. GratsHEr: (i) Report of the Comm. on Mathematical 
Tables, 1873, p. 31-33, 150; (ii) ‘“‘On circulating decimals with special reference to Henry 
Good ’s ‘Table of circles’ and ‘Tabular series of decimal quotients’ (London 1818-1823),” 
Cambridge Phil. Soc., Proc., v. 3, 1878, p. 185-206; (iii) “On a property of —— frac- 
tions,” Phil. Mag., s. 5, v. 7, 1879, p. 321-336. We have sivntie fe) 
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Statement of 1861 [MTAC, v. 2, p. 87] concerning Goodwyn: “His manuscripts, an 
enormous mass of similar calculations, came into the ion of Dr. OLintHus GrEGory, 
and were purchased by the Royal Society at the sale of his [Gregory's] books in 1842.” 
In the above publications of Glaisher it is stated that no trace of the papers could be found 
at the Royal Society. Neville’s statement, p. xi, “Goodwyn left a mass of papers, no one 
knows what became of them,” is therefore slightly misleading. 


882[F].—BencT Persson, “On a Diophantine equation in two unknowns,” 
Arkiv Mat., v. 1, 1949, p. 45-57. 


The table on p. 57 gives all solutions of the diophantine equation 
x +x+ (p+ 1)/4 = y* for all values of the prime p < 100 of the form 
4n + 3, except the value » = 3 which is handled separately on p. 45 of 
the paper. 
Ivan NIVEN 
University of Oregon 
Eugene, Oregon 


883[F].—B. vAN DER Pot & P. Spezia, “The primes in k(p),”” K. Nederl. 
Akad. Wetensch., Proc., v. 54 A, 1951, p. 9-15, 1 plate. 


The primes of the quadratic field k(e***/*) are presented in graphical 
form. Each integer in the field is assigned to the center of a hexagon. The 
totality of integers thus fills the plane. Those integers which are primes are 
indicated by black hexagons, the non-primes by white hexagons. The plate 
inserted between p. 13 and 14 gives a graphical representation of all primes 
in the field whose norms are less than 10000. The beautiful symmetrical 
diagram contains 7416 black hexagons out of approximately 36276 hexag- 
onal cells giving a density of primes of about one out of five. 

An enlargement of the central part is given in another figure. The paper 
contains a table of solutions (c, 8) of the quadratic partition 


p= b6m+1=2+ 3f for p < 10000. 
There are two misprints, communicated by the author: 


for p = 1559 read 1549 
for p = 7459 read 7549 


This table was used to construct the diagram. The authors were un are 
that previous tables of this extent exist. Jacosi’s table' for p < 1’ 47 is 
more than a century old. Other tables by REUSCHLE? and CuNNI am? 
extend to p = 13669 and p = 125683 respectively. [See also UN ‘4. ] 
A comparison of Cunningham’s table with the present one r <als no 
discrepancy. 

A similar diagram for the primes of the field k(z), based on squares, was 
published by van der Pol.* 

D. H. L. 


1K. G. J. Jacost, “Uber die Kreistheilung und ihre Anwendung auf die Zahlentheorie,” 
Jn. reine angew. Math., v. 30, 1846, p. 166-182. 

2K. G. REuSCHLE, Mathematische Abhandlung, enthaltend: Neue Zahlentheoretische 
Tabellen, Stuttgart, 1856. 

3A. J. C. CunnincHam, Quadratic Partitions. London, 1904; Quadratic and Linear 
Tables. London, 1927. 

4B. vAN DER Pot, Verslagen van de Maatschappij Diligentia, The Hague, 1946. A copy 
of this diagram woven in red and white squares hangs on the wall of the reviewer's study, 
a gift from the author. 


140 RECENT MATHEMATICAL TABLES 


884[H].—W. P. Tuiysen, ‘Over de geitereerde machtsverheffing,” Simon 
Stevin, v. 27, p. 177-192, 1950. 


The table on p. 185 gives approximate solutions of the equations: 
(1) t = x! 
(2) t = x 


for given x. Equation (1) has two solutions:!¢ = A,u;1<A<e<yu< om. 
Similarly equation (2) has two solutions: t = 6,7; 0<a<e'<7< 1. 
5D values (mostly) are given for \, and 4D values for o, r, all for x = 0(.01).06, 
e~* = .06599, .1(.1)1, 1.2, e/* = 1.447. The values of u corresponding to 
x = 1, 1.2, e'/* are ~, 14.77, e respectively. 

On p. 184 there are graphs! illustrating A, u, o, 7 as functions of x. 


R. C. ARCHIBALD 
Brown University 


Providence, R. I. 

1 The curve ¢ = x‘ has been frequently studied before. For example: (a) J. F. C. HEssEL, 
“Uber das merkwiirdige Beispiel einer zum Theil punctirt gebildeten Curve das der Glei- 
chung entspricht y = Nx," Archiv Math. Phys., v. 14, 1850, p. 169-187; H. ScHEFFER, 


“Uber die ~~ die Gleichung y = Nx dargestellten Curven,” Archiv Math. Phys., v. 16, 
1851, p. 133-137. 

the role which ¢ = w plays in Cantor’s theory of transfinite numbers will be recalled; 
see, for example, G. Cantor, Math. Annalen, v. 49, 1897, p. 242-246 (also the English 
translation of P. E. B. Jourpatn, Open Court, 1915, p. 195-201). 


885(I].—H. E. Sawzer, “Coefficients for polar complex interpolation,” 
Jn. Math. Phys., v. 29, 1950, p. 96-104. 


This article is supplementary to two earlier papers,! both of which are 
devoted to the computation of the coefficients of interpolation formulas for 
the representation of a function f(z) defined at m successive points 2, where 
k ranges from [$n — }] to [42], the points being spaced equally along the 
arc of a circle about the origin of coordinates. 

In the first of these papers the interpolation formula was written in 


the form 
f(z) ~ XL Li (P) f(z), 


where the point Pywas defined to be P = (2 — 20)/(z; — 20) = p + ig, and 
where P" was represented by , + ign. Explicit formulas were given for 
the real and imaginary parts of the coefficients L,, which are functions 
not only of p, and gp, but also of 6, since z = pe. The three cases of 3-point, 
4-point, and 5-point interpolation were treated. 

In the second paper the interpolation formula appeared as follows: 


f(z) ~ (X arf) /X a, 


where we write 


Q,= (ze — 2) *** (Se — Se-1)(Se — Sear) *** (Se — Sn) (2 — z) 


If new quantities A,“ = (z — 2,)a, are introduced, and if 2 = 0, % = 1, 
= 1+ e*, 23 = 1+ ec” + ce?” 2, = —e-# — ¢-2, etc., then A,™ can 
be computed as functions of @. This article gives the explicit values for 
these functions for 3-point to 9-point interpolation inclusive. 
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The article under review provides tables of the functions L;™ and A,™, 
computed in the first case to 8D and in the second case to 9D, with some 
doubt as to the ninth place, over the range @ = 1°, 5°, 10°, 15°, 20°, and 30°. 

H. T. Davis 


Northwestern Univ. 
Evanston, IIl. 

1H. E. Sauzer, “Formulas for direct and inverse interpolation of a complex function 
tabulated along equidistant circular arcs,” Jn. Math. Physics, v. 24, 1945, p. 141-143, 


MTAC, v. 2, p. 73, and “‘Alternative formulas for direct interpolation of a complex function 
tabulated along equidistant circular arcs,” Jn. Math. Physics, v. 26, 1947, p. 56-61. 


886(K].—R. L. ANDERSON & T. W. ANDERSON, “Distribution of the circular 
serial correlation coefficient for residuals from a fitted Fourier series,”’ 
Annals Math. Stat., v. 21, 1950, p. 59-81. 


The model for observations x;, i = 1, 2, ---, N, here used is 
Xi — wi = p(X — win) — Hi, 


with x9 = x, and po = ua, in which the {u;} are random disturbances each 
independently and normally distributed with zero mean and variance ¢?, 
and 4; is a linear form in Fourier terms. The coefficients of the Fourier 
terms are the usual regressions of the x; on these terms. The fitted values 
of the y; are denoted by m; and the serial correlation coefficient 


N N 
R= D¥ (x — mi) (xin — ma) / & (x: — m,)* 
i=1 i=l 
(with lag 1), where mo = m,, is used to test hypotheses concerning p. 
Under the null hypothesis of zero serial correlation (p = 0), Table I, p. 65, 
gives to 3D the exact 5% and 1% points for R for P periods used in the fitted 
series for the following cases: P = 2, N = 6(2)60: P = 2,4, N = 8(4)100; 
P = 2, 3, 6, N = 12(6)150; P = 2, 12/5, 3, 4, 6, 12, N = 24(12)300. 
Table II, p. 66-67, gives exact 5% and 1% significance points for R (the 
distribution of R is asymmetric in these cases) for single periods for P = 3, 
N = 6(6)150; P = 6, N = 12(6)150; P = 4, N = 8(4)100, 108(12)144; 
P = 12, N = 12(12)300. 
. ¢.<. 


887(K].—D. G. CHAPMAN, “Some two sample tests,” Annals Math. Stat., 
v. 21, 1950, p. 601-606. 


Let s represent the difference of two independent variables, each of 
which has a #-distribution with » — 1 degrees of freedom. Table I contains 
values of P(O = s S so) ; the values of m and so considered are : m = 2(2)12, ©; 
So = .5(.5)8(2)12, 21, 30, 50, 100. Table II contains values of so such that 
P(|s| 2 so) = .01, .05 for m = 2(2)12, «. 

Another place where 1 percent and 5 percent points for the absolute 
difference of two independent /-statistics can be obtained is in the 1948 
edition of FisHzerR & YATES’ tables.! Let ¢, have a ¢-distribution with 
degrees of freedom while ¢, is independent of ¢; and has a ¢-distribution 
with m2 degrees of freedom. Let d = (t; — t2) cos 6. Table V; (by P. V. 
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SUKHATME) contains values of dy) such that P(|d| > do) = .01, .05 for 
mM, M, = 6, 8, 12, 24, and @ = 0°(15°)90°. 


J. E. Wats 
U. S. Bureau of the Census 


Washington, D. C. 


1R. A. FisHer & F. Yates, Statistical Tables for Biological, Agricultural, and Medical 
Research. London, 1948. 


888[K].— JEROME CORNFIELD & NATHAN MANTEL, “Some new aspects of 
the application of maximum likelihood to the calculation of the dosage 
response curve,”’ Amer. Stat. Assn., Jn., v. 45, 1950, p. 181-210. 

Table 1 (p. 185-188) gives Z/P, —Z/Q, Z?/Q?—ZY/Q, and ZY/P+2Z?/P?, 

all to 5D, for X = 5(.01)10, and for X = 0(.01)5 by symmetry. X — 5 = Y 

is a normal deviate, Z being the corresponding ordinate; P is the proba- 

bility of a deviate less than Y, and Q that of a larger deviate. The purpose 

of the table is to minimize the arithmetic of the maximum likelihood esti- 

mate of the probit regression. There are detailed instructions for calculation. 
H. W. Norton 


University of Illinois 
Urbana, Illinois 


889[I, K].—D. B. DeLury, Values and Integrals of the Orthogonal Poly- 
nomials up ton = 26. Toronto, Published for Ontario Research Founda- 
tion by University of Toronto Press. 1950, v + 33 p. 17.5 & 25.2 cm. 
$1.25. 


The polynomials tabulated are those of CHEBYSHEV' and Gram and 
are the FISHER & YATES? polynomials denoted by £,’(x) = &,’(x, m). 

The values of é,’(x,) are tabulated for r = 0(1)n —1 and for 
n = 3(1)26. 

The sums of squares of the é’s are given at the foot of each column. 

For r < 14 and m < 26 values of 


i, = = &,' (x, n) dx 
0 


and 


I= fee, nde 
=f 


are tabulated. 

The tables presented in the book were prepared for the purpose of 
furnishing an arithmetical basis for a numerical integration procedure 
appropriate to situations in which the observed or measured ordinates are 
subject to a random error. 

Of considerable interest is the discussion concerning the fitting of a 
function in two variables by means of orthogonal polynomials. 

Two examples of the use of the tables are given. The first deals with 
evaluating an integral over a range of equally spaced data. The second 
example shows how the total number of trees on a given area can be esti- 
mated from counts on square plots situated centrally on square blocks 
within the area. The orthogonal polynomials are used to determine the 
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equation of the surface. The equation is then integrated to yield the desired 
estimate. 

It is of interest to point out that a table of orthogonal polynomials 
and their derivatives has been prepared by W. F. Brown, Jr. & C. W. 
Dempsey [MTAC, v. 4, p. 224-5]. 

A comparison between the values of the polynomials tabulated by 
DeLury and by Brown & Dempsey show no discrepancies. It is to be 
regretted that the values of the ratio of ¢’ to = are not given in the book. 

The book is clearly printed. It should be a welcome addition to the 
library of statisticians and mathematicians. 


Jack SHERMAN 
Beacon Laboratories 


Beacon, N. Y. 


Further values of these orthogonal polynomials have been given by 
ANDERSON & HousEMAN? for m = 3(1)104, r = 1(1)min(5, ”) and by VAN 
DER REYDEN‘ for m = 5(1)13, r = 1(1) various values less than 9 and for 

= 14(1)52, r = 1(1)9. 

The introduction describes the uses of the table in fitting polynomials 
by least squares, in finding integrals of the resulting polynomials, in carrying 
out the associated analyses of variances and in determining confidence 
intervals in one and in two independent variables. 

Values of &,’ for r < 5 are said to be reprinted from FisHer & YATEs.? 
The author states that ‘The director of the Onderstepoort Veterinary 
Research Institute and Dr. D. vAN DER REYDEN*‘ also have kindly granted 
permission to reproduce portions of the tables which appeared in the 
Onderstepoort Journal.”” The present tables are, however, correct in two 
places where minus signs are missing from van der Reyden’s table (for 

= 17, r=5, &' = —3 and m = 18, r=5, &’ = —1. Vander Reyden 
gives values of the polynomials for x < #). And a common factor in the 
values for = 25, r = 8 in van der Reyden’s table has been removed. 

For more complete discussions of these orthogonal polynomials see the 
review by E1rsENHART [MTAC, v. 1, p. 148-150] and papers by Brrce,' 
WEINBERG® and van der Reyden‘ and the references given in these places. 
For punched card procedures using such tables as the one under review 
see a paper by LILA KNUDSEN.’ 

K. J. ARNOLD 


University of Wisconsin 
Madison, Wisconsin 


1P. L. CHEBysHeEv, “Sur I’interpolation par la méthode des moindres carrés,” Akad. 
Nauk, Leningrad, Mémoires, s. 7, v. 1, no. 15, 1859, p. 1-24. Oeuvres, v. 1, p. 471-498. 
2R. A. FisHer & F. YATEs, Statistical Tables for Biological, A gricultural and Medical 
Research. Third ed., London, 1948. (First ed., 1938.) 
3? R. L. ANDERSON & E. E. HoUSsEMAN, Tables of Orthogonal Polynomial Values Extended 
to N = 104. Ames, Iowa, Iowa State College of Agriculture and Mechanic Arts, 1942. 
(Research Bulletin 297, p. Arey Agricultural Experiment Station, Statistical Sections.) 
LMTAC, v. 1, p. 148-150, v. 814 
4D. VAN DER REYDEN, ‘ ‘Coles Glee by the orthogonal polynomials of least squares,” 
Onder stepoort Journal of Veterinary Science and Animal Industry, v. 28, 1943, p. 355404. 
5 RayMonpD T. BirGE, ‘‘Least-squares’ fitting of data by means of polynomials,” Review 
of Modern ~ a v. 19, 1947, p. 298-347. 
EINBERG, “Mathematical appendix,” ibid., p. 348-360. 
ma F. KNUDSEN, “A punched card technique to obtain coefficients of orthogonal 
polynomials Am. Stat. Assn., Jn., v. 37, 1942, p. 496-506 
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890[K].—A. K. Gayen, ‘‘The distribution of the variance ratio in random 
samples of any size drawn from non-normal universes,” Biometrika, 
v. 37, 1950, p. 236-255. 


The author is interested in the non-normal sampling distribution of the 
two statistics: w = the ratio of two estimates of the variance obtained in a 
one-way classification for the analysis of variance (homogeneity), v = the 
variance ratio of two independent samples (compatibility). To introduce 
non-normal populations the EDGEWORTH series, 


fle) = o(e) —*2 (x) +  go(e) + gC), 


is used. With an Edgeworth series population the author evolves the 
frequency densities of the two statistics. These distributions are then 
related to their corresponding ‘‘normal”’ sampling distributions and correc- 
tion terms for the skewness and excess are determined. The coefficients of 
As (= VB) and Ay (= 6: — 3) are tabulated in Table 4 (p. 252-255) to 4D 
for the 5 percent points of Fisher’s Z for the degrees of freedom », = 1(1)6, 
8,12, 24, ©; m = 1(1)6, 8, 12, 20, 24, 30, 40, 60, 120, «o. These tables enable 
one to compute the error one would make in using the ordinary tables 
associated with the normal population distribution of Z for the 5 percent 
level if he were aware of the A; and \, of his population and was willing to 
assume an Edgeworth type distribution. In this connection the author 
notes that the Edgeworth restriction would not be serious if the size of 
sample was sufficient to enable one to neglect terms of order N-*. If \3 and 
\4 are unknown one might use their sample estimate in the formulas to get 
an indication of the size of error involved. 


C. F. Kossack 
Purdue University 


Lafayette, Ind. 


891[K].—N. F. GJEDDEBAEK, “Contribution to the study of grouped obser- 
vations. Application of the method of maximum likelihood in case of 
normally distributed observations,” Skandinavisk A ktuarietidskrift, v. 32, 
1949, p. 135-159. 


Let 


o(2) = (2x)-Fexp (— 2/2), @@) =f" ody, o@) = 


o(x + y) — (x) _ (x+y) — ¢'(x) 
(x+y) — 8x)’ *~ Se + y) — Be) * 


Table 1, p. 151-154, contains 4D values of 2:(x,y) for y = 0(.1)4 and 
x = —.1h(.1)4, where h is the largest integer < 5y + 1. Table 2, p. 155- 
158, contains 4D values of 22(x, y) for the same values of y and x in Table 1. 
Given y, both z, and z, are symmetrical about the value x = —y/2. Table 3, 
p- 159, contains 4D values of 2; and 2 for y+ and x = —4.5(.1)0. The 
functions 2; and z: are used in an iterative procedure for determining the 
maximum likelihood estimates of the mean and standard deviation of a 


dey(x) 
dz ° 








a(x,y) = — 
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normal population when the observations from this population have been 


grouped. 


Bureau of the Census 
Washington, D. C. 


J. E. Wats 


892(K].—A. HaLp, The Decomposition of a Series of Observations Com- 
posed of a Trend, a Periodic Movement and a Stochastic Variable. 
Copenhagen, G. E. C. Gads, 1948, 134 p., 17.2 & 25.1 cm. 8 kr. 


The author defines polynomials £,(x,;), of degree r=1, ---,1(v=1,-+-,m; . 
i = 1,---,k), by means of the relations £,.. = 0 (a dot in a subscript indi- 
cates a mean with respect to the corresponding variable), » = 1, ---, J, and 


n k 
LL (kroi — Eri) (Ewe — Eni) = Oforu <1,7 = 2,+++,);4 = 1,-+-,) —1. 


y=l i=l 
Here £,,; is used to denote £,(x,;): 

The leading coefficients are assigned the value unity and then the remaining 
coefficients are determined by the defining equations. For k = 1 these 
become the ordinary CHEBYSHEV polynomials over equally spaced abscissas. 
The tables give values generally to 6S, of these polynomials of degrees 
1 to 5 for the combinations k = 3, m = 5,7; k = 5, nm = 5,7, 8, 9, 10, 11. 
For k = 4, m = 25 the degrees run from 1 to 4 and for k = 12, n = 10, 
the degrees run from 1 to 6. In addition explicit formulas for the poly- 
nomials tabulated are given (the author was unable to derive the general 
explicit formula) and the values, again generally to 6S, of the elements of 
the three matrices ||£,.;||, ||S,,||, ||.S’-.|], for the polynomials tabulated in which 


n k 
Srp ” 2 L (Erni - £505) (Envi - fi), (Sr =0,7r# ), 


and 


n k 
Sp = LL bribes. 


y=l i=l 


¢c.6. 
893[K].—H. O. Hart Ley, ‘“‘The maximum F-ratio as a short-cut test for 
heterogeneity of variance,’’ Biometrika, v. 37, 1950, p. 308-312. 


The ratio of the largest to the smallest in a set of k variance estimates 
is proposed as a short-cut test in place of BARTLETT’s test,' 


M = N log. {N- 9 ms} - z Vt log. sé 
t t 


where the estimate s? is based on » degrees of freedom and N = >> ». 
t 


Since log, s? is approximately normal with variance 2/(v — 1), 
Frax(a) = exp {we(a)V2/(» — 1)} 


gives approximately the 100a% point of the ratio when the & estimates 
are based on the same number, », of d.f. and where w(a) is the corre- 
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sponding point of the distribution of the range in samples from normal 
(PEARSON & HARTLEY’). Tables 1 and 2 give comparisons of the approxi- 
mate and exact 5% points of Finax to 3S fork = 2andv = 2(1)10, 12, 15, 20, 
and for vy = 2 and k = 2(1)12, respectively. The two sets of exact results 
were used to obtain second approximations by the adjustment, Fyrex(a) 
(second approx.) = Fmax(a) (first approx.) (1 + qq), with g, and q fitted 
to, the exact values. The second approximations, with exact values in the 
first row and the first column, are given in Table 3, which shows upper 
5% points of Fnsx to 3S in a set of k mean squares all based on » d.f. for 
k = 2(1)12 and vy = 2(1)10, 12, 15, 20, 30, 60, «. 
Leo Katz 

Michigan State College 
East Lansing, Mich. 

1M. S. Bart ett, “Properties of sufficiency and statistical tests,’ R. Soc. London, 
Proc., v. 160A, 1937, p. 268-282. 


2E. S. Pearson & H. O. Hart ey, “The probability integral of the range in samples 
of m observations from a normal population,” Biometrika, v. 32, 1942, p. 301-310. 


894[K].—H. O. Hartiey & E. S. Pearson, ‘‘Tables of the x?-integral and 
of the cumulative Poisson distribution,’ Biometrika, v. 37, 1950, p. 
313-325. 


Let P(x?, v) = zine} f ex held x, 
x2 


The authors have tabulated (p. 318-325) P(x?, v) to 5 places of decimals 
for v = 1(1)20(2)70 degrees of freedom and x? = 0(.001).01(.01).1(.1)2(.2)- 
10(.5)20(1)40(2)134. At the same time the table provides the values of the 


c—1 
cumulative Poisson distribution since P(x?, v) = }> e~™m*/i! with m = x?/2, 
i=0 
c = v/2. Hence for m = 15 the complete Poisson distribution is provided 
and for m > 15, only the truncated Poisson sum up to c = 35 (v = 70). 
Methods of interpolation in the table are thoroughly discussed. The authors 
point out that Pearson’s table of the incomplete I'-function’ and MoLtna’s 
table of the Poisson exponential limit? are more extensive, the former for 
x? and the latter for the Poisson distribution. But Pearson’s table involves 
a troublesome transformation and Molina’s tables do not give P(x?, v) for 
odd values of v. In addition it may be noted that Salvosa’s table of the 
Type III function, both areas, ordinates and the first six derivatives? may 
be used to evaluate P(x?, v) by setting ¢ in Satvosa’s table = (x? — v)(2v)}, 
and a;? = 8/v. A comparison by the reviewer with the tables of P(x?, 2) 
given in the Elements of Statistics by Davis & NE son‘ for small values 
of x? in a few instances revealed no divergencies, but for v = 2, x? = 5 
Davis and Nelson give .08208 49986, while the present table gives .08209. 
The same phenomenon was noted for v = 9, x? = 3.0. No doubt the table 
will be found useful in statistical investigations. 
L. A. AROIAN 
Hughes Aircraft Company 
Culver City, Calif. 


1 KarRL PEARSON, Tables of the Incomplete T'-Function. London, 1922. 
2E. C. Mouina, Poisson's Exponential Binomial Limit. New York, 1945. 
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?L. R. Satvosa, ‘Tables ¢ pre's type III function,” Annals Math. Stat., v. 1, 
1930, Ee 191-198. Appendix, p. 1 
. T. Davis & a Elements of Statistics with Applications to Economic 
pn ‘acioleaiend Ind., 1935. 


895[K].—J. M. Howe t, “Errata to ‘Control chart for largest and smallest 
values,’ ’’ Annals Math. Stat., v. 21, 1950, p. 615-616. 


In the paper cited in the title [MTAC, v. 4, p. 207-208] the author 
notes that numerical errors exist in Table I (as was pointed out by MTAC’s 
reviewer). Using values appearing in a paper by Gopwin,! HOWELL now 
gives a corrected version of Table I to 4D. 

<<. 


H. J. Gopwin, “Some low moments of order statistics,” Annals Math. Stat., v. 20, 
1949, p. 279-285. [See MTAC, v. 4, p. 20. 


896[K].—R. F. Linx, “The sampling distribution of the ratio of two ranges 
from independent samples,” Annals Math. Stat., v. 21, 1950, p. 112-116. 


Consideration is given to a sample of m ordered observations drawn from 
a population with variance o”. Let w = (x, — x,)/o. Then for two samples 
of sizes m, and m, drawn from populations of the same variance, the ratio 
w;/W2 becomes the ratio of the two ranges. Table II gives values of R for 
all combinations of m, and mz = 10 and for a = .005, .01, .025, .05, and .10, 
such that 
Pr (w,/wz < R) = 


where the samples are drawn from populations with the same variance. 
These values may be used for testing the hypothesis that two independent 
samples were drawn from normal populations with the same variance. This 
test is therefore comparable to the F test. 


T. E. BICKERSTAFF 
University of Mississippi 
University, Miss. 


897[(K].—E. Lorp, ‘‘Power of the modified /-test (u-test) based on range,” 
Biometrika, v. 37, 1950, p. 64-77. 


For a sample of m independent observations from a normal population 
with unknown standard deviation, o, the u-test is the ratio of the normal 
deviate to an independent range estimate of the standard error. The power 
of a two-tailed test, 8(p), against an alternative p = £/c, is the sum of the 
two power components, §’(p) = Pr {lower tail|p} and 8’’(p) = Pr {upper 
tail|p}. The standardized error of NEYMAN & TOKARSKA! is pz such that 
1 — B(p.) = @ for any significance level a. Six tables are given to compare 
the standardized errors of the u-test and the ¢-test fora = .05 and a = .01. 

Table I, p. 67, shows to 5D the power components 6’(p) for p = 1(1)10 
and #’(p) for p = 1(1)3 and » = 1(1)20, © of the ¢-test, where » is the 
number of degrees of freedom for estimating variance. Table II, p. 68, 
gives the standardized errors for the ¢-test, p.os to 3D and p. to 2D for 
v = 1(1)20, 30, 60, 120, ~. 

Table III, p. 70, gives the power components for the u-test to 5D for 

= 2(1)20, ©; B’(p) for p = 1(1)10 and @’(p) for p = 1(1)3. Table IV, 
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p. 71, shows the standardized errors for the w-test, p.o5 to 3D and p.o, to 2D 
for m = 2(1)20, o. 

In Table V, p. 72, there are given values of p.os to 2D for the u-test 
based on the average range in m samples of independent observations 
each, for m = 2(1)20 and m = 1(1)10(5)20, 30, 60, 120. From Table V, it 
appears that the optimum grouping of N (= mn) observations is one such 
that 6 = m = 10 in general. 

Table VI, p. 76, gives the ratio of the standardized errors at 5% and 1% 
levels of the u-test to the #-test for 3D for m = 2(1)20, ©. The ratio is 
= 1.024 for the 5% level and =1.045+ for the 1% level. 

Leo Katz 
Michigan State College 
East Lansing, Mich. 


1J. NevMaNn & B. ToxarskA, “Errors of the second kind in testing ‘Student’s’ hy- 
pothesis,” Am. Stat. Assn. Jn., v. 31, 1936, p. 318-326. 


898[K].—P. B. Patnaik, ‘The use of mean range as an estimator of vari- 
ance in statistical tests,” Biometrika, v. 37, 1950, p. 78-87. 


Let x1, X2, °**,%n be a random sample of m observations, arranged in 
ascending order of magnitude, drawn from a normal population with mean u 
and variance o”. Define the range by w, = x, — x;, and if there are m such 
independent samples with m observations each, the mean range by Wm, a. 
Further d, is defined by E(w,) = od,. The author approximates the dis- 
tribution of ,,, by means of the first two moments of ,,, and the x? 
distribution with v degrees of freedom and scale factor c. In Table 1, needed 
for the distribution of #,,,, are given the values of d, for m = 3(1)10 to 
4D, v (the degrees of freedom which are not integers) to 3D, c, the scale 
factor to 4D, for m = 1(1)5, and m = 3(1)10. Table 2 gives the compari- 
sons of the exact and approximate percentage points of u = d,x/Wm,, for 
n = 5,8,10, m = 1,3,5,10, a = .1, .05, .02, .01, where x is distributed 
normally with mean zero, and variance o? and a denotes the percentage 
points. The exact results are from Lorp’s computations,' while the approxi- 
mate results are based on the author’s approximation to the u distribution. 
Table 3 compares the exact power function of u with the author’s approxi- 
mate Student’s non-central ¢, for w/o = 1, 2,3, u #0, m = 5,10, m = 3,5 
for a one-tailed region, a = .05, to 3D. For the symmetrical or double- 
tailed test, comparisons of the power are made for m = 1, n = 10; m = 2, 
n=10;m=4,n=5;m=5, n = 3, and m = 2, nm = 5 usually to 3D. 
The results are particularly useful in quality control and wherever it is 
desirable to reduce computation time by use of the mean range in place of 
the sample standard deviation as an estimator of the population standard 
deviation. 

L. A. AROIAN 
Hughes Aircraft Company 
Culver City, Calif. 
1 E. Lorp, “The use of range in place of standard deviation in the ¢-test,’’ Biometrika, 


v. 34, 1947, p. 41-67; ‘Power of the modified ¢-test (u-test) based on the range,” Biometrika, 
v. 37, 1950, p. 64-77. [RMT 897] 
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899[K].—L. W. Po.tak, assisted by U. N. Ecan, All Term Guide for 
Harmonic Analysis and Synthesis Using 3 to 24; 26, 28, 30, 34, 36, 42, 
44, 46, 52, 60, 68, 76, 84 and 92 Equidistant Values. (Department of 
Industry and Commerce, Meteorological Service, Geophysical Publica- 
tions, v. 2.) Dublin, Stationery Office, 1949, xx, 185 p. 24.5 & 30.2 cm. 
£2 s2. 


The purpose of the present work is to facilitate the work of fitting the 
Fourier series 


y = Po t+ picosx + p2.cos 2x + --- + qsinx + q@sin 2x + --- 
to a set of m equidistant values, where m is any one of the numbers given 
above in the title. 

The tables are supplementary to those published by the author in 1947 
and described in MTAC, v. 2, p. 306. The tables of the “Schedule,” which 
comprises most of the work (185 p.), give the values to 5D of the functions 

2xkm . 2aekm 


cos ' sin , 
n n 








where m is any one of the numbers given above in the title, m < n, and 
k < 4m. 

In a typical entry 2 is printed in large blackface type, let us say » = 20, 
and for this number there are then 10 tables corresponding to k = 1 to 
10 = 4n. The value of & is printed in smaller blackface type at the top of 
each table. Two arguments are given, one for m, denoted by i by the author, 
which ranges over the positive integers, and the other for i’ = km — sn, 
where s is chosen so that km — sn is positive and less than n. 

Two other entries are included, which refer to tables in the author’s 
Rechentafeln zur harmonischen Analyse, Leipzig, 1926, by means of which 
the multiplication of the harmonic terms is facilitated. The choice of m given 
above was guided by the choice of values found in the Rechentafeln. 

The present work is a valuable addition to other similar tables which 
have been computed and which are listed in MTAC, v. 1, p. 193. 

H. T. Davis 


Northwestern Univ. 
Evanston, Ill. 


900[K].—E. M. Scurock, Quality Control and Statistical Methods. New 
York, Reinhold Publishing Corporation, 1950. xi + 213 p., 17.1 XK 24.1 
cm. $5.00. 


Tables 9.7a, p. 69-71, and 9.7b, p. 73-122, of this book were prepared 
with the apparent purpose of routinizing the construction of control charts 
for industrial applications covering most situations which are likely to be 
encountered in practice, in so far as problems dealing with the per cent 
defective are concerned. 

Table 9.7a labeled, ‘“Table of Width of 3-Sigma Band in Per Cent,” 
gives the amounts to be added to and subtracted from the average per cent 
defective or process level, #, in order to obtain 3-sigma upper and lower 
control limits for sample per cent defectives. Thus Table 9.7a gives values 
of 3[p(1 — p)/N] to 4D, where N is the sample size, for N = 10(10)100- 
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(100)4000(200)5000(500) 10000, and # = .001, .0025, .005, .01, .015, .02(.01)- 
-1(.02).2(.05).5. For non-tabulated j’s and N’s, the author states, “although 
the relationships are curvilinear, straight line interpolation between rows 
and columns may be used for all practical purposes.” 

Table 9.7b, Section One, gives 3-sigma upper and lower control limits 
for the number of defectives to the nearest integer and for the per cent 
defective to 2D and also gives the expected number of defectives for ranges 
of # and N as above. GRANT! gives a comparable table of upper and lower 
control limits alone. Section II of this table is a rearrangement of the same 
information for. numbers of defectives, giving separately for each N = 100- 
(100)4000(200)5000(500)10000, the values for the same j’s as above up to 
.35. Thus Section One of Table 9.7b with different sample sizes on each 
page is convenient for testing control with respect to an established or 
standard process level of per cent defectives, whereas Section II of Table 
9.7b “‘is extended for use in situations where several different quality char- 
acteristics, each with a different quality level, will be checked on samples 
of the same size.”’ 

Tables 9.7a and 9.7b occupy some 54 pages of the book and because of 
the wide ranges of sample sizes and values of per cent defective should be 
of considerable value as a ready reference to the quality control engineer 
for industrial problems involving per cent defective. Figures 9.7a and 9.7b 
are included in the book for obtaining limits for sample sizes other than 
those given in Tables 9.7a and 9.7b or for values of per cent defective larger 
than 35%. These figures can also be used for 3-sigma confidence intervals 
for the per cent defective of the lot or material samples. Examples are given 
for the use of Tables 9.7a, 9.7b and Figures 9.7a, 9.7b. 

F. E. GRUBBS 


Ballistic Research Laboratories 
Aberdeen Proving Ground, Md. 


1E. L. Grant, Statistical Quality Control. New York, 1946. 


901[K].—C. A. Wi.iams, Jr., “On the choice of the number and width 
of classes for the chi-square test of goodness of fit,” Am. Stat. Assn. 
Jn., v. 45, 1950, p. 77-86. 


This paper contains a presentation and discussion, intended mainly for 
the non-mathematical user of tests of goodness of fit, of results of a paper 
by MANN and WALD! which make it possible to optimize, in a certain sense, 
the number of classes and the choice of class-boundaries for the chi-square 
test. Tables are given showing, for the 1% and the 5% significance level 
and the sample sizes N = 200(50)1000(100)1500, 2000, the optimum num- 
ber & of classes, as well as some information related to the power of the test. 
The author suggests that the tabulated values of k may be halved with 
little loss of power. 


Z. W. BIRNBAUM 
University of Washington 


Seattle, Wash. 


1H. B. Mann & A. WALD, “On the choice of the number of class intervals in the appli- 
cation of the chi-square test,” Annals Math. Stat., v. 13, 1942, p. 306-317. 
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902[L].—AKaApDEmMIIA Nauk, SSSR. Matematicheskif Institut im. V. A. 
Steklova, Institut Mekhaniki, Tabliisy znachenit Funkisii Besselia ot 
mnimogo Argumenta. Edited by I. M. VinoGrapov & N. G. CHETAEV, 
Moscow, Leningrad, 1950, xii, 404 p. 39 roubles. 


Here we have finally the work which was started in 1930 and ready for 
publication in 1941 [MTAC, v. 3, p. 66 and v. 5, p. 112], after collaboration 
of various individuals and academies, as set forth in the preface. The world 
war delayed publication, but after further revision by the present editors, 
an edition of 2,500 copies was finally printed last year. 

After the introductory matter there are two sections (p. 1-201, 203-403) 
each containing four tables. These tables are of (1) Io(x) to 8D; (2) Ih(x) to 
9D; (3) 24-*Ko(x) and (4) 27'K,(x) to 8S or 8-12D; (5) I;(x) and (6) I_4(x) 
to 8S; (7) Ko(x) and (8) K,(x) to 8S or 8-12D; all for x = 0(.001)10, A. 

The only publications referred to in the literature list are: Watson, 
A Treatise on the Theory of Bessel Functions, 2nd ed., 1944 and the Russian 
translation of 1949. Gray-MATHEWsS—MACRoBERT, A Treatise on Bessel 
Functions. The date given is 1931, but presumably 1936 was intended. 
A Guide to Tables of Bessel Functions as in MTAC, v. 1, no. 7, 1944. FMR, 
Index, 1946. And finally there is a reference to tables of Jo(x) and J;(x) in 
BAAS, Report, 1893 and 1896, but these tables are attributed to ALpts, 
and not to the real author A. LopGE; these 9D tables were for x = 0(.001)- 
5.1, A. This mention of the name Aldis, where there was no reason for it, 
and general considerations, suggest very definitely that other unmentioned 
important tables were well known to the compilers and computers. Never- 
theless, it is undoubtedly true that most of the values given in the volume 
before us are new. 

Let us recall just what was available before, apart from the tables of 
Lodge. In R. Soc. London, Proc., v. 64, 1899, p. 218-223, ALpis gave tables 
of In(x), h(x), for x = [.1(.1)6; 21D], [6(1)11; 18D]. BAASMTC, Math. 
Tables, v. 6, 1937, contained tables of Jo(x), Ii(x), for x = [0(.001)4; 8D], 
[4(.001)5; 7D], 6? and of Ko(x), Ki(x), for x = [0(.01)2; 8D], [2(.01)4; 9D], 
[4(.01)5; 10D], 6%. In the case of Jo, J; it was found that the Lodge tables 
were not dependable in final digits, and hence the calculations were derived 
by subtabulation in the 21D values of Aldis. In NBSMTP, Table of the 
Bessel Functions Jo(z) and J,(z) for Complex Arguments, 1943, second edition 
1947, there are tables of Ip(x), J:(x), for x = [0(.01)10; 10D]. These tables 
occur again in NBSCL, Table of the Bessel Functions Y,(z) and Y,(z) for 
Complex Arguments, New York, 1950, published almost simultaneously with 
the volume under review. In this volume for x = [0(.01)10; 10D] are also 
tables of (3) and (4) 2x-'Ko(x), 2x-'Ki(x). The only previous tables of 
Ix;(x) are those of A. N. Dmynix, 1915 [MTAC, v. 1, p. 287], I(x), for 
x = [0(.1)15;4-6S]; I_4(x), for x = [0(.1)7;5-6S]. But Dinnik’s other 
tables have been found notoriously unreliable [MTAC, v. 2, p. 379]. 
Finally there is NBSAMS 1, Tables of the Bessel Functions, Yo(x), Y:(x), 
Ko(x), Ki(x), 1948 [MTAC, v. 3, p. 187-188], Ko(x) and K(x) for x = [0- 
(.0001).033(.001)1; 7S], A*. 

From these statements the exact relation of the Russian volume to what 
had appeared before (which excludes NBSCL, 1950) is clear. The only 
exact duplication is in Ip(x), [i(x), for x = 0(.001)4 in the BAASMTC 
volume, and in the Lodge J;(x) table. 
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Thus, if the tables in this new volume are found to be accurate, 
notable addition has been made to the large amount of tabular material 
dealing with Bessel functions. 

NOTE ADDED IN PROOF: A friend (who desires that his name be withheld) 
has drawn my attention to ten additional errata, each of 3 to 10000 units, 
in final decimal places: 


Function x for read 
In(x) 1.009 175 075 
2.46 8582 7582 
T(x) 4.06 6 3 
4.67 6 3 
4.87 4 1 
4.88 5 2 
5.07 3 0 
5.08 6 3 
To(x) 6.79 9981 1981 
9.73 46664 56664 


Because of these, and many other less serious errors already observed these 
Russian tables of the Academy of Science must be characterized as decidedly 
unreliable. 


R. C. ARCHIBALD 
Brown University 
Providence, R. I. 


1 Qn spot-checking the 9D table of J;(x) in the work under review with the 10D table 
of NBSMTP, 1943, one difference of 4 units in the end-figure was found at x = .5, and 
49 differences of one or two units in end-figures were found for the following values of x: 
OG, Ss, Bi, BP, .1, ts, B, 1.82, 1.03, 1.11, 1.15, 1.16, 1.17, 1.18, 1.25, 1.3, 1.35, 1.36, 
1.37, 1.39, 1.5, 1.55, 1.57, 1.59, 1.7, 1.85, 1.87, 1.94, 1.97, 1.98, 2.67, 2.68, 2.69, 2.7, 2.85, 
3.03, = 3.95, 4.35, 4.55, 4.65, 4.71, 4.72, 4.76, 4.77, 4.79, 4.85, 5.05. 

At ouary ones of these 50 places, except for x = .5 (because of the very weakly printed 
left part of Lodge’s ending 4, it was copied in the Russian table as a 1), the values of the 
Russian table were identical with the values of the Lodge (1893) 9D table of J:(x). Hence, 
it seems evident that the Lodge table was copied completely without any checking what- 
ever. Since the British Committee has told us that the 9D tables of I(x) and J,(x) ‘“‘were 
subject to errors of two or three units in the last decimal, so that rounding off would not 
yield 8-figure values of the standard of accuracy set for the Committee’s tables,” we can 
assert that the Russian tables are probably incorrect at each of the places indicated above. 
Furthermore, spot-checking the 8D table of Jo(x), differences were found at x = .07, 3.3, 
4.4, where there was evident rounding off of Lodge values. Since correct 8 D values are given 
in the BAASMTC tables it is definitely ——— that the Russians also copied their table 
of Io(x), up to x = 5.1, rounded off from Lodge, and that they did not have the Jo(x) table 
in BAASMTC. Since the Lodge table of Tie) abridged to interval .01 was reprinted in 
Gray & MATHEWS, i, 1922, - all the. apebelite errors noted above are also to be 
found in G. & M. and G. M. & M 

In a brief spot checking met y x= 5.1 up to x = 7, unit differences in end figures 
were found at x = 5.3, 5.45, 5.57, 6.15, 6.35, 6.45, 6.75. 


903[L].—P. A. Carrus & C. G. TREUENFELS, “Tables of roots of incom- 
plete integrals of associated Legendre functions of fractional orders,” 
Jn. Math. Phys., v. 29, 1951, p. 292-299. 


P,}(cos @) is the associated Legendre function of order one and fractional 
index yu. Tabulated are the zeros of 


d 
P,3(cos @) = 0, 76 P,(cos 6) = 0 


treated as a function of the index. 








fu 
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Table I gives the first 50 values of u for @ = 90°(5°)175°, while Table II 
gives the first 50 values of »v for @ = 90°(5°)130°. The &-th root is denoted 
by ux, %. In both tables entries are mostly to 5S. Values were derived from 
power series, for the smaller roots, and from asymptotic expansions due to 
Pat! for the larger roots. Convergence difficulties prevented extension of 
Table II beyond 130°. 

No indication as to accuracy of the tables is given, though it is stated that 
interpolation of the roots 6-wise is practicable with five point Lagrangean 
interpolation. The reviewer has spot checked the tables by differencing and 
in a number of cases found that third differences are not smooth, indicating 
that some errors are present. For example, in Table I, entry for k = 1, 
@ = 130° is undoubtedly in error. The entry 1.3001 should probably read 
1.2975. In Table II, there is a typographical error for the entry k = 49, 
@ = 90°: for 88.000 read 98.000. In Table II, for & large, third differences 
are not smooth. The reviewer has not differenced all the entries. In absence 
of error, interpolation with four point interpolation should be sufficient to 
insure full accuracy of Tables I and II. 

In Tables III and IV, values of the integral 24 Pn'(x) }2dx are pre- 
sented. Here x = cos 0, m stands for yw or v, and tabulations are given for 
each value of the index found in Tables I and II. However, in Table III, 
6 = 90°(5°)145°, in Table IV, @ = 90°(5°)130°. Entries are to 4S and 5S. 
Finally, Table V gives values of fo'{P,(x)}*dx to 5D for m = 0(0.05)1.0 
together with second differences. 

All tables are new, though Pat' had previously determined some yu and » 
values of the functions considered for various orders and values of @ in the 
first quadrant. 


YupDELL L. LuKE 
Midwest Research Institute 


Kansas City, Missouri 

1B. Pat, “On the numerical calculation of the roots of the equations P,"(«) = 0 and 
z P,™(u) = 0 regarded as equations in m,”” Calcutta Math. Soc., Bull., v. 9, 1918, p. 85-95 
and v. 10, 1919, p. 187-194. 


904[L].—A. H. Carter & O. A. WILLIAMs, JR., “A new expansion for the 
velocity potential of a piston source,’ Acoustical Soc. Amer. Jn., v. 23, 
1951, p. 179-184. 

A small table is given (p. 182) of the values of the hypergeometric 
function F(—(m — 1), —m;1;6*) which the authors denote by fo». These 
functions are actually polynomials of degree m—1 in 6. Thus 


fu=1+20, fo = 1+ 6b? + 304. 
The 6D tables give fon(b) for 
m = 1(1)8 and b = 0(.1)1. 


There is also a table of the coefficients of fom for m = 1(1)8. These 
functions occur as factors of coefficients in the expansion of the velocity 
potential of a piston source in terms of HANKEL functions. 


D. H. L. 
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905[L, P].—R. A. CLark, “On the theory of thin elastic toroidal shells,”’ 
Jn. Math. Physics, v. 29, 1950, p. 146-178. 


Table I, p. 156, gives 4D values of 


% {sin #/(1 + A sin #)}#dt 


for \ = 0(.1).5 and @ = 0°(2°)90°. 
Table II, p. 167, gives 3D values of the real and imaginary parts of 

T(x) and dT(x)/dx for x = 0(.05)5 where T(x) is that solution of 

"ta 

yi xT = 1 
which is asymptotic to 7/x for large (positive) x. In terms of LOMMEL’s 
function 

T(x) = — $(ix)*Soa[3(tx) 4). 
A. E. 


906[L].—R. A. FisHErR, “Gene frequencies in a cline determined by selection 
and diffusion,” Biometrics, v. 6, 1950, p. 353-361. 


This paper contains a table of the solution g(x) of the non-linear differ- 
ential equation g’’ = 4xq(1 — q) subject to the two point boundary condition 
q(0) = 3, g(~) = 0. The table gives g(x) to 8D for x = 0(.02)5.7. Beyond 
this range g(x) < 10-*. The table was computed on the EDSAC [MTAC, 
v. 4, p. 61-65]. 

D. H. L. 


907(L].—CaARL-ERIK FROBERG, ‘On determination of proton-proton inter- 
action from scattering experiments,” Arkiv for Fysik, v. 3, 1951, p. 1-23. 


The CovuLoMB wave equation containing an additional nuclear inter- 
action potential V(r) is 


(1) d*y/dr? + [k? — 2akr- — l(l + 1)r7] y = — Vir)y. 


If V(r) = 0 a system of two linearly independent solutions can be given by u%, 
and x, such that ~,(0) = 0, u, ~ sin (kr + ) asr > and » ~ cos (kr + 7) 
as r—, with 7 = —alog (2kr) — 3/9 + arg (ia +/1+ 1). If V(r) is different 
from zero an additional phase-shift 5(%) appears in the asymptotic expression 
for the regular solution y of (1), so that y(0) = Oand y ~ sin [kr + 9 + 4(k) ] 
as r >. The author discusses in this paper the determination of V(r) if 
5(k) has been determined experimentally for every positive k. The relation 
between V and 4 is given by 

(2) V(r) = - : (2ak)? f a~* sin 5(a) ¢(a, A) da 

where 


ola») = (£ us(o)oa) 


p=A/ (2a) 


and 


X = 2akr. 
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To determine V(r), a knowledge of the function g(a, A) is necessary. 
Therefore the author gives asymptotic expansions for u and v. The main 
part of the paper consists of two numerical tables. The first table is for 
%;(p) and v'(p) with / = 0 for a = .01 and O < p < 6, a = .02, .03 and 
0<p <4, a= .04 and 0 < p < 2, a = .06(.02).12(.04).2(.05).3(.1)1 and 
0 < p < 1: the intervals in p vary. The second table is for g(a, A), a ranging 
from .01 to © and A from .012 to .120, with various intervals. The last 
section of the paper gives the computation of V(r) for two cases, one where 
sin 6 = const. a? and a = .02, the other where sin 6 is approximately con- 
stant (.79) and .006 < a < .02. Three smaller numerical tables are included. 
The tables were computed on the Swedish automatic computing machine 
BARK in Stockholm [MTAC, v. 5, p. 29-34]. 


Maria A. WEBER 
California Institute of Technology 


Pasadena 4, California 


908[L].—R. GramMEL, “Tafeln der verallgemeinerten Kreisfunktionen 
Sin(4)v, Cos(4)v, Sin(6)v, Cos(6)v, Sin(8)v, Cos(8)v,”’ Ing.-Arch., v. 18, 
1950, p. 251-254. 


The functions x = cos (m)v and y = sin (m)v have been defined in an 
earlier paper! as the inversion of the integrals 


(1) o= f fall) ea f fall) - 
where 
f.() = (1 — t*)9o-!*, 


For » = 2 these functions reduce to trigonometric functions. For any even 
integer m, they provide a parametric representation of the curve 


(2) x+y = 1, 
and they have certain periodic and symmetry properties with half-period 


T, = f fa(t) dt. 


In the present paper 6D numerical tables of cos (m)v and sin (m)v are 
given for nm = 4,6,8 and v = 0(.01)1(.1)4 and for 20/7, = 0(.05)1. The 
earlier paper’ gives 3D values of x, for m = 2(1)10, ©, and the present 
paper 6D values of $7, for m = 4, 6,8. For m > 2 one has x < m, < 4. 

The infinite series of the earlier paper! have been used for the computa- 
tion of the functions for 0 < v < 4,, symmetry and periodic properties 
combined with interpolation, for v > }2,. The original computations were 
carried to 7 or 8D, and the author states that all six decimals printed are 
correct. 

Differencing of the tables was used as a check, the relation (2) was 
checked for every fifth entry, and for m = 4 the connection with Jacobian 
elliptic functions was used to check the tables. 

The computations were carried out by F. JINDRA. 

A. E. 


1R. GRAMMEL, “Eine Verallgemeinerung der Kreis- und Hyperbelfunktionen.” Jng.- 
Arch., v. 16, 1948, p. 188-200. 
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909[L].—E. Kreyszic, ‘Uber den allgemeinen Integralsinus Si(z, «),”’ Acta 
Math., v. 85, 1951, p. 117-181. 


The function 
(1) Si (2, u) = f t-* sin t dt 
0 


was introduced by A. WALTHER.! It is related to the incomplete gamma 
function in the same manner as the ordinary sine integral, Si (z, 1), is related 
to the exponential integral function. Together with the generalized cosine 
integral, 


(2) Ci (2, uw) = f t-* cos ¢ dt, 
0 


the integrals investigated in this paper account for the incomplete gamma 
function, sine-, cosine-, and exponential integrals, and Fresnel’s integrals. 

The author investigates Si and Ci in the complex domain. z = x + iy, 
p=a+ isp, 0 <a <2 in (1), and 0 < a < 1 in (2). In the theoretical 
part he discusses power series expansions in ascending powers of 2, the 
connection with other functions including confluent hypergeometric func- 
tions, asymptotic expansions as z+, and a convergent expansion in 
incomplete gamma functions. He investigates the zeros and the product 
representation of the functions in question and comments on the numerical 
computation of the zeros. 

Table 1 (p. 155-6) gives 3D values of Si (x, a) for x = 0(.2)4(.5)20, 
@ = .25(.25)1.75. 

Table 2 (p. 157-163) gives 2D (sometimes 3S) values of the real and 
imaginary parts of Si(x + iy, a) forx = 0(1)20,y = 0(1)5, a = .25(.25)1.75. 

Table 3 (p. 164) gives 2D values for the real and imaginary parts of the 
first three zeros of Si (z, a) for a = 0(.25)1.75. 

Table 4 (p. 165) gives 5D values of Si (, a) fora = .05(.05)1.95, 1.99, 
1.995, 1.999, 2, and 5D values of Ci (©, a) for a = .05(.05)1. 

Table 5 (p. 166) gives 8D values of 


Pla + 1,9) = [ewar 


for z = 1, —1, 4 (in the last case real and imaginary parts) and m = 0(1)10. 
There are relief diagrams and charts illustrating the functions tabulated, 
other diagrams; and an extensive bibliography of papers on both theory 
and applications, and also of numerical tables of the sine integral and 
FRESNEL’s integral. 
The tables are said to be accurate to within a unit of the last decimal. 
Although the relation 


Ci (s, uw) = z* sins — w Si (2,4 + 1), 0<a<i, 


allows one to compute Ci from Si, tables of Ci (z, 4) would be a useful 
addition. Together with the tables of the present paper they could serve to 
compute the incomplete gamma function in the complex domain. 

A. E. 


1 A. WALTHER, “‘Anschauliches zur Gibbsschen Erscheinung und zur Annaherung durch 
arithmetische Mittel,”” Math. Zeit., v. 42, 1937, p. 355-364. 
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910[L].—-L. M. MILNE-Tuomson, Jacobian Elliptic Function Tables. Dover 
Publications Inc., 1950. xii, 132 p., 12.7 K 18.9 cm. $2.45. 


The purpose of this book is explained in the opening sentence of its 
preface. ‘The widespread belief that calculations involving elliptic func- 
tions are difficult is due not to the nature of the calculations themselves 
but to the lack of suitable numerical tables wherewith to perform them.” 
In fact, a glance at A. FLETCHER’s ‘Guide to Tables of Elliptic Functions” 
[MTAC, v. 3, p. 229-281 ] suffices to show that while there is a considerable 
body of numerical material for the use of a professional computer, and also 
for some special purposes, yet there is no handy and easily available book 
of tables which could be used by the general practitioner of the various 
sciences as easily, and as safely as, say, a table of logarithms. The present 
book attempts to fill this gap. 

The standard notations are used, except that the parameter m = k? is 
used rather than the modulus k, sn (u|m) is written for sn (u,k), and 
similarly for cn, dn. The complementary parameter is m, = 1 — m. 

The first 39 pages are devoted to the definition of Jacobian elliptic 
functions and a description of their principal properties, examples showing 
the use of the tables, complete elliptic integrals, Fourier series and power 
series expansions of elliptic functions, periods, zeros, poles and residues, 
formulas for special values of the argument, change of argument (including 
double and half arguments and addition theorems), change of parameter 
(transformations), approximations, complex arguments, differentiation, 
WErersTRASS’s g-function, integrals of Jacobian elliptic functions, elliptic 
integrals and the zeta function, conformal mapping, factorisation of cubic 
and quartic polynomials, and application to the pendulum problem. In 
brief, these 39 pages give practically all the information constantly needed 
when using elliptic functions. 

The following tables and graphs form the main body of the book under 
review. 

sn (u|m). Graphs for m = 0, 3, 1 and 0 < u < 2x (p. 41). 5D tables 
with first differences for m = 0(.1).5 and w = 0(.01)2, m = .6, .7, .8 and 
u = 0(.01)2.5, m = .9 and uw = 0(.01)3, and m = 1 and u = 0(.01)3(.1)6.5, 
with values of the quarter-period K at the foot of each column (p. 42-61). 

cn(u|m). Graph (p. 63) and tables (p. 64-83) in the same range as for 
sn, except that for technical reasons the values for m = 1 (which are common 
to cn and dn) and u > 3.5 are distributed over various pages of the dn-table. 

dn (u|m). Graph (p. 85) and tables (p. 86-105) again for the same values 
of m and u as for sn, except that dn(u|0) = 1 for all « and need not be 
given, and the column thus saved is used to accommodate tables for 
cn(10u!1) = dn(10u|1) up to the point (w = 1.29) where these functions 
vanish (to 5D accuracy). 

Complete elliptic integrals and the nome (p. 106-109). 7D tables of 
K, K’, E, E’, and 8D tables of g = exp(—7K’/K), and q, = exp(—7K/K’) 
for m = 0(.01).5, with values of m, printed on the right margin and an 
arrangement like that of trigonometric tables to enable the user to read off 
values for .5 < m < 1. 

Z(u). 7D table (p. 112-123) for m = 0(.1).6 and u = 0(.01)2, m = .7, 8 
and « = 0(.01)2.5, and m = .9, 1 and uw = 0(.01)3, all with A”. 
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The book is very well printed, and the legibility of the tables is excellent, 
except that the table for Z(u) seems to have been reproduced photograph- 
ically from the Proceedings of the Royal Society of Edinburgh (v. 52, 1931, 
p. 239-250) and the considerable reduction in size makes for small figures 
and a somewhat crowded page. 

A. E. 


911[L, V].—M. Munk & G. RAWLING, Tables of Chaplygin Functions. Naval 
Ordnance Laboratory Memorandum 10816, 57 pp. (1950). [Available 
to government contractors only.] 


This paper consists of three tables of certain functions designated by 
the generic title of CHAPLYGIN functions, which arise as particular solutions 
of Chaplygin’s differential equation. They are closely related to the physical 
plane stream function or to the LEGENDRE reciprocal potential of the hodo- 
graph plane in the computation of plane two-dimensional compressible 
flows. The adiabatic index is taken as 1.4. 

This differential equation is: 

2 
r(r — 1) 222 + Cm + 1 + Be — (m + 1)]S* — $6m(m — 1)Zm = 0 





in which 6 = (y — 1)! = 5/2, and 7+ is a dimensionless speed variable 
related to the Mach number by r = M*(28 + M7)“ and to the reduced 
velocity w by r = w’. 

Particular solutions tabulated are Z, = F(adm, bm; Cm, 7), Where 
Am + bm = m + B, Ombm = —4$Bm(m — 1), Cm = m+ 1; and Z_n, where 
Z_m denotes a second independent solution; mot negative values of the 
parameter m. This second solution Z_,, contains a logarithmic term and 
has been normalized to approach unity as w tends to zero. 

The other functions were computed from the relations 


Ln = W"Zm L_n = W"Z_m, 
Vn = p(mLn — mwL,,), Vem = p(mL_m — mw) 


where p = (1 — w*)5/?, and a prime denotes d/dw. 

Values of the functions for w = 6—}, sonic speed, were found by inter- 
polation using NEwTon’s backward formula. 

The tables give 6S values of the functions 


Bans pte ee ie ) ee | ol | Ries Wm) Ym 


for m = 2(1)10 and w = .01(.01).5 and for w = 6~! to which latter value 
of w a separate table is devoted. 

The tables were produced by the computational facilities of the Naval 
Ordnance Laboratory. 


G. A. RAWLING 
U. S. Naval Ordnance Laboratory 


Silver Spring, Md. 


912(L].—NBSCL, Table of the Bessel Functions Yo(z) and Y,(z) for Complex 
Arguments. Columbia Univ. Press, 1950. xl, 427 p. 20 X 26.5 cm. $7.50. 
This volume maintains the high standards of the earlier tables, pre- 
pared by the computation group under A. N. Lowan. It presents values 
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of the real and imaginary parts of the second solutions Yo and Y; of 
the Bessel equation, according to the magnitude p of the argument pe‘* 
over the range [0(0.01)10; 10D] and according to the phase angle » 
over the range [0°(5°)90°]. Auxiliary tables of VY, — (2/2)Jy log p and 
Yi — (2/x)Ji log p + (2/xp)e-**, over the range [p = 0(0.01)0.50; ¢ = 0°- 
(5°)90°; 10D], are provided, to facilitate interpolation for small values of p. 
In addition, the position of some of the zeros of Yo = Yj’ and of Y; and Y,’ 
off the real axis are given, together with values of Y, and Y; at some of 
these points. A table of 5-point Lagrange interpolation coefficients is also 
provided, for use with the tables. In the introduction are given contour 
plots for the general behavior of Yp and Y; in the complex plane; and a 
discussion of the properties of the functions and the method of calculation 
of the tables. 

Together with the earlier volume,’ giving corresponding values of Jy 
and J; for complex values of the argument, these tables will allow many 
acoustical and microwave problems to be computed for complex boundary 
conditions on cylindrical surfaces. It is particularly satisfactory to have 
both functions J of the first kind and functions Y of the second kind, for 
the same values of the order and argument. Many applications of Bessel 
functions require knowledge of both functions and so many tables of these 
functions, heretofore published, have lost a great part of their value because 
only the J’s were tabulated, not the Y’s. In the present case the computing 
group is to be congratulated in resisting the temptation to compute J2, J; 
etc., and in turning to the more difficult but more useful Y, and Yj. 

Perhaps some of our larger digital machines will eventually take the 
place of tables, making up in speed of output what they inevitably lack in 
judgment; they show little sign as yet of their ability to produce further 
tables, as carefully: planned, edited, and produced as the tables under 
review. Perhaps, in the distant future, when we have a large digital com- 
puter in each laboratory, we will not need tables of the sort reviewed here 
(though the present reviewer thinks we always will). At any rate, for a 
long time to come, these tables and the many others produced by Lowan’s 
able group will be standard and much used tools of physicists and engineers. 


P. M. Morse 
Mass. Inst. of Tech. 


Cambridge, Mass. 
1NBSCL, Table of the Bessel Functions Jo(z) and J,(z) for Complex Arguments. 2nd 
ed., New York, 1947 [MTAC, v. 3, p. 25]. 


913(L].—A. T. Price, ‘Electromagnetic induction in a semi-infinite con- 
ductor with a plane boundary,” Quart. Jn. Mech. Appl. Math., v. 3, 
1950, p. 385-410. 


Table I, p. 398, gives 3 and 4D values of Yi, Ye, Z:, Z: for —, 7 = 0, .5, 
1(1)5, 10 where 


Y, = f fae cos un du, Z2.™= f fre sin un du, 
0 0 


fi=i1- uv2{(1 + u*)? — u2}}, fo = 2u? — uv2{(1 + u*)t + ut}. 
A. E. 
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914[L].—P. Ruopes, ‘‘Fermi-Dirac functions of integral order,”” R. Soc. 
London, Proc., v. 204A, 1950, p. 396-405. 


The functions 


Fala) = [xe + 1) 


arise in the theoretical treatment of assemblies of particles subject to 
Fermi-Dr1rac statistics. For the cases » = 3, } numerical tables are known.! 
In the present paper ” is a positive integer. A series convergent for 
negative 7, a relation between F,(n) and F,(—7), and a polynomial approxi- 
mation for large (positive) 7 are derived. 
Table 1 (p. 404) gives the approximating polynomial for m = 1(1)4. 
Table 2 (p. 404) gives 7D values of F,,(n)/n! for m = 1(1)4, 7 = —4(.1)0. 


A. E. 


1J. McDoucaLt & E. C. Stoner, “The computation of Fermi-Dirac functions,” 
R. Soc. London, Phil. Trans., v. 237A, 1938, p. 67-104. 

C. TRUESDELL, “On a function which occurs in the theory of the structure of polymers,” 
Annals Math., s. 2, v. 46, 1945, p. 180 [MTAC, v. 1, p. 445]. 


MATHEMATICAL TABLES—ERRATA 


In this issue references have been made to Errata in RMT 883 (VAN DER 
Pot & SPEZIAL!), 889 (DELURY), 894 (HARTLEY & PEARSON), 895 (HOWELL), 
902 (AKADEMIA Nauk SSSR). 


190.—A. J. C. CUNNINGHAM, Binomial Factorisations. V. 2, London, 
1924. 


Page 189, line 6 
for 39249421 read 30249421 = 1291-23431 
D.H.L. 


191.—A. J. C. CUNNINGHAM & H. J. WooDALL, Factorisation of (y" + 1). 
London, 1925. 


Page 17, » = 66, delete the factor 3 
n = 77, insert the factor 463. 


N. G. W. H. BEEGER 
Nicolaas Witsenkade 10 
Amsterdam C, Holland 


192.—J. P. STANLEY & M. V. WILKEs, Table of the Reciprocal of the Gamma 
Function for Complex Argument. Computation Centre, Univ. of Toronto, 
1950 [MTAC, v. 5, p. 25-26]. 


The following errors occurred in the preparation for press. A number 
of illegible entries are also noted. Thanks are due to Miss C. M. Munford 
of the University Mathematical Laboratory, Cambridge, and to Dr. van 
Wijngaarden of the Mathematical Centre of Amsterdam, who helped in the 
discovery of these errors. 
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Argument 
U x ¥ Column* 
—0.17 0.18 read 179387 
—0.39 0.23 read 330349 
—0.38 0.31 read 380853 
+0.43 0.38 for 510159 read 519159 
+0.35 40-.44 (arg) for +.035 read +-0.35 
+0.02 0.44 1 read 84522 
—0.41 0.45 for 525223 read 535223 
+0.43 0.45 read 535157 
+0.33 0.45 for 485888 read 385888 
+0.02 0.64 2 for 210157 read 201157 
+0.46 0.72 2 for 706845 read 706849 
+0.43 0.79 for 597722 read 697722 
—0.01 84.88 (arg) for 0.01 read —0.01 
—0.45 0.85 for 1—.441265 read — 1.441265 
0.00 0.92 1 Insert “‘ —0” 
+0.04 0.97 for 393497 read 393947 
Argument 
V x y Column* 
—0.19 0.06 read 42836 
+0.33 0.06 for 70200 read 70290 
+0.42 0.16 1 for 185696 read 185686 
—0.45 .20-.24 (arg) for —.045 read —0.45 
—0.35 .20-.24 (arg) for —0.30 read —0.35 
-0.15  .32-.36 (arg) for 0.15 read —0.15 
+0.50 0.45 read “ +-” 
+0.10 0.46 read “ +-"" 
+0.05 0.47 read “‘ +" 
+0.26 0.52 1 for 605490 read 695490 
—0.40 -56-.60 (arg) for +0.40 read —0.40 
—0.31 0.64 2 for 435506 read 433506 
—0.09 0.71 for 844477 read 884477 
—0.40 .72-.76 (arg) for 0.40 read —0.40 
—0.35  .72-.76 (arg) for 0.35 read —0.35 
—0.25 .72-.76 (arg) for 0.25 read —0.25 
—0.45 0.80 2 for +1.278102 read +-0.278102 
+0.34 0.88 1 for 1422458 read 1442458 
—0.33 0.90 read 807496 
—0.34 0.92 2 read 820268 
—0.15 0.93 read 1352942 
—0.27 0.95 for 1155247 read 1115247 
—0.31 0.99 for 1012513 read 1102513 
* Some columns are repeated on two pages, “1” identifies the first repetition, ‘‘2"’ the 


second. 
C. G. GOTLIEB 
Computation Centre 
Univ. of Toronto 
Toronto, Ontario 
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UNPUBLISHED MATHEMATICAL TABLES 


122[B].—H. E. Sauzer, Table of n!/x"*!. Manuscript in possession of the 
author, NBSCL. 


This table gives, to within a unit in the 9th place, the values of the 
function n!/x*+! for nm = 0(1)10 and x = .1(.1)10. The table was computed 
with the aid of PETERs & STEIN’s values of the powers of the reciprocals 
on p. 36-57 of the Anhang to J. PETERS, Zehnstellige Logarithmentafel, v. 1. 


123[F].—A. GLopEN, Factorization of 2n?+ 1 for m = 1(1)1000. Type- 
written manuscript of 11 leaves deposited in UMT FILE and Scientific 
Computing Service, London, and in possession of the author, 11 rue 
Jean Jaurés, Luxembourg. 


A table extending to = 3000 is contemplated by the author. 


124[F].—A. GLopEN, Tables of Quadratic Partitions p = a? + 3b’ for 
100000 < » < 200000. Typewritten manuscript of 18 leaves, deposited 
as in UMT 123. 


The argument ? is a prime of the form 6m + 1. The table gives a and b 
for each p with some 35 omissions which are promised in an addendum 
[for previous tables up to = 125000 see RMT 883 ]. 


125(F].—R. M. Rosinson, Tables of Integral Solutions of |y* — x*| <x 
for x < 10°/9. Tabulated from punched cards and deposited in the UMT 
Fite. Cards in possession of the author, University of California, 
Berkeley. 


There are 5 tables. The first gives all solutions arranged according to x, 
of which there are 1242. Table II lists the 332 cases in which y’ = x’. 
Tables III is arranged in order of |y* — x*|. Table IV is for y? > x* and 
Table V is for y* < x*. Tables IV and V omit the trivial solutions x = 4k? + 1, 
y = 8k® + 3k, which are starred in Tables I and III. 


126[L].—G. Firaxe & Y. L. LuKeE, Tables of the Function (1 — uZ)-'*. 
Lithographed manuscript, 4 leaves, deposited in UMT FILE and avail- 
able from Midwest Research Institute, 4049 Pennsylvania Ave., Kansas 
City, Mo. 

The table gives 8D values of this function for 


z= .7(.01)1.25, « = .09(.01).12. 


127(L].—C. P. Green, J. H. Litre & D. W. Rayno.ps, Extensive Tables 
of the Exponential Integral. A thesis, Chemical Engineering Dept., Univ. 
of Tennessee. Knoxville, 1951, xxv + 111 leaves. 


The main table is of —Ei(—x) for x = 15(.0001)16. There are also 
tables of Ei(x) and —Ei(—x) for x = 15(1)50 and interpolation coefficients 
as introduced by Coutson & Duncanson [Phil. Mag., s. 7, v. 33, 1942, 
p. 745-760]. 
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128[L].—G. Jones & D. Urrorp, Table of the Functions C(P) = K,(P)/ 
[Ki(P) + Ko(P)] and of PC’(P). Lithographed manuscript, 4 leaves, 
available asin UMT 126. 


K, and Ko are the usual Bessel functions. The tables give C(P) and 
—PC'(P) to 7D for P = 0(.002).1(.01).3(.02)1(.1)2(.5)10(10)100, and also 
for P = .35(.05).95. 


129[L].—Y. L. Luxe & D. Urrorn, Tables of the Function Ko(x) = o*Ko(t) dt. 
Lithographed manuscript, 3 leaves, deposited as in UMT 126. 


The table gives 8D values of Ko(x) and of the auxiliary functions A,(x) 
and A,(x) defined by 


Ko(x) = (In 2 — y — In x)A;(x) + A2(x) 
for x = 0(.01).5(.05)1. 


130(L].—UNIvERsITY OF TORONTO COMPUTATION CENTRE, Tables of Spheri- 


cal Bessel Functions for Semi-imaginary Argument. Photo copy, 2 leaves 
deposited in UMT Fine. 


The tables give 8S values of the real and imaginary parts, absolute 
values and arguments of 


(2x/a)—te-*t4 I na(xe™?) 
(2x/a)—te-**"4 Vn44(xe**?) 
for nm = 0,1, 2,3; x = 0(1)10. 


AUTOMATIC COMPUTING MACHINERY 


Edited by the Staff of the Machine Development Laboratory of the National Bureau 
of Standards. Correspondence regarding the Section should be directed to Dr. E. W. 
Cannon, 415 South Building, National Bureau of Standards, Washington 25, D. C. 


Notes on Numerical Analysis—5 
Table-Making for Large Arguments. The Exponential Integral 


The evaluation of a function defined by a definite integral, for the 
complete range of argument — © to +, is usually performed in several 
stages. For small and moderate values of the argument x the integral is 
evaluated by means of an ascending series in powers of x, or perhaps by 
numerical quadrature. For very large values of x, numerical values are 
obtained by means of an asymptotic series. 

The exponential integral, for example, defined by the equations 


Ei(x) = f * let dt 
(1) 
~Ei(—x) -f lett 
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has been tabulated by the use of ascending series, 


n 





Ei(x) = y+ log2+ 5 
n=1 


n-n! 
(2) i 
—Ei(—x) = —y—-—logx + Db =< 


n=l 


n-n! 
for a range of x up to about 15, and from the asymptotic series 


(3) Ei(x) ~~a-te*(1 + 1b! + 2x? + ---) 
—Ei(—x)~x-te*(1 — 1!x-1 4+ 2b? — ---) 


for larger values of x. 

In some cases solutions for moderate values of x are obtained by numeri- 
cal solution of the differential equation satisfied by the function, though this 
method does not appear to have been used for the exponential integral. 

The use of the ascending series is rather cumbersome for |x| > 10, and 
the asymptotic series may not, without further refinement by the ‘‘con- 
verging factor’’ method of ArrEy,' provide the required number of figures. 
This note gives a method for extending the tabulation, without any use of 
asymptotic series, to cover the whole range of x from —© to +. The 
method is not restricted to the exponential integral, but the latter provides 
a useful illustration. 

The extension to large values of x is most easily performed by using 
an argument z = 1/x, and it is easy to show that the auxiliary function T, 
defined by the relation 


(4) T = €*Ei(x) — x7 = e""Ei(z) — 2 
satisfies the differential equation 
(5) T” — (4 -—22°)T+ 2° = 


where dashes denote differentiations with respect to z = 1/x. 
We can also show that the Taylor series at z = 0 has the form 


(6) T= 2+ 2+ --- + mig 4+ --- 


and is divergent for all z > 0. In spite of this fact we can use finite-difference 
equations and numerical integration to solve equation (5). 

We assume that values of T are available, from previous computation 
with the ascending series, at z = + 0.1 (x = + 10), and we shall fill in 
the values for z = —0.1(.0i1) + 0.1 by the use of relaxation methods, 
producing by this means a table from which both Ei(x) and —Ei(—x) can 
be quickly obtained from the relation 


(7) Ei(i/z) = e/#(T + 2) 


For numerical purposes equation (5) is replaced by a set of difference 
equations to be satisfied at every pivotal point in the range of integration 
and typified by 


Ti + Ty — {2 + W(eq-* — 2a0-*)} To + Ae? + A(To) = 0 
A(To) = (—64/12 + 68/90 — ---) T) 


(8) 
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At z = 0, T is zero and no further equation is needed, though it is easy 

to show from (6) and (8) that the final solution should satisfy the equation 
Ti+ Ti - 2h? + A(T>) = Oats = 0. 


The set of equations (8) is ideal for relaxation or any iterative process, 
as can be seen from the numerical values of the coefficients of 7, for 
h = 0.01, which are given in the following table: 


2 coeff. of T> zZ coeff. of T> 

01 9802 —.01 10202 

02 602 — .02 652 

03 118 + 4/81 — .03 132 + 70/81 

.04 37 + 15/16 — .04 44 + 3/16 

05 16 + 2/5 —.05 19 + 3/5 

.06 8 + 64/81 — .06 10 + 52/81 

.07 5 + 1397/2401 — .07 6 + 1796/2401 
.08 4 + 13/256 — .08 4 + 213/256 
.09 3 + 1639/6561 —.09 3 + 5239/6561 


The relaxation process follows familiar lines. A first approximation T® 
is obtained by neglecting A, an approximation to which is computed from 
the differences of T. The “‘missing’”’ differences near the ends of the range 
are filled in by the methods of the previous note. (See MTAC, v. 5, p. 92-95.) 

Two applications of the iterative process, the first of which produced a 
maximum change in T of 200, the second of 15 in the ninth decimal place, 
gave the results shown in the table appended. 

From these values we can obtain by subtabulation a table from which 
Ei(x) and Ei(—x) can be computed for any |x| > 10. The only other table 
known to exist in this range is that of Coutson & DUNCANsON,? in which 
interpolation is by no means trivial. 

It is also interesting to consider other numerical methods. The powerful 
method VII of Fox & Goopwin,' for example, replaces equations (8) by 


(1 — h?f,/12)T, + (4 — 2 f_4/12)T_1 — (2 + Sh*fo/6)To 
(9) + h?(Sgo/6 + gi/12 + g-s/12) + A(To) = 0 
A (To) = 6°/240 — 1368/15120 + --- 
where 
f =z‘ — 2s, g=s". 


Only two cycles of the iterative process are now required ; but the coeffi- 
cients of 7; and 7_,; are no longer unity, and heavier arithmetic is called 
for in the relaxation. Equations (9) have been applied as a check, however, 
and indicate that the last figure given is nowhere in error by more than 
one unit. For small values of z, furthermore, more significant figures could 
be obtained trivially, a result of the large coefficient of 7) in the finite- 
difference equations. 

We also considered the solution by step-by-step methods of the first 
order equation 
(10) T’ —2°T+1=0. 


Method II of Fox & Goodwin® produces the recurrence relation 


(1 — $hay*)T, = (1 + $ha)T —h+A=0 


(11) A = (—6*/12 + 85/120 — 87/840 + ---)Tj. 
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Starting with T = 0 at s = 0 and integrating outwards, we run into trouble 
associated with the fact that the coefficient of 7; passes through zero. 
Other methods have similar disadvantages, and indeed the presence of the 
large coefficients of T) in equations (8) suggests immediately that step-by- 
step methods will be difficult, relaxation relatively easy. 

If we apply step-by-step methods using Taylor series, we soon run into 
trouble in working outwards from z = 0, owing to rapid growth of the 
complementary function introduced by rounding-off errors. Since it is, 
however, possible to work inwards from z = 0.1, this method was used as a 
spot check. 

An independent computation has been performed by T. VicKERs, who 
used asymptotic series and interpolation in existing tables to produce values 
of the function 

S = xe*Ei(x) = xT + 1 


in our previous notation. The limitations of the asymptotic series make it 
difficult to obtain more than seven decimals in S, effectively one less than 
we obtained in T. Vickers points out, however, that the differences of S 
converge more rapidly than those of T; in particular linear interpolation 
in S can be performed with an error of less than 2 units in the fourth 
decimal place. 

Tables of S and T, with facilities for interpolation, are given in Table 1. 

L. Fox & J. C. P. MILLER 


NBSCL 
TABLE 1 
st 10° T Smt aa S bn2 
—.10 8436666 123964 0 0.9156333 1992 
—.09 6936157 129007 0 0.9229316 2107 
— .08 5564725 134430 0 0.9304409 2238 
—.07 4327803 140290 0 0.9381742 2379 
— .06 3231260 146634 1 0.9461457 2533 
—.05 2281454 153538 1 0.9543709 2711 
—.04 1485301 161068 1 0.9628675 2904 
— .03 850351 169333 1 0.9716550 3124 
—.02 384890 178446 1 0.9807555 3375 
—.01 98058 188554 1 0.9901942 3663 
.00 0 199858 1 1.0000000 3994 
01 102063 212593 2 1.0102063 4383 
.02 417046 227108 2 1.0208523 4848 
.03 959551 243875 3 1.0319850 5414 
.04 1746477 263623 4 1.0436619 6124 
.05 2797795 287587 7 1.0559559 7080 
.06 4137929 318250 10 1.0689655 8467 
.07 5798138 358767 11 1.0828305 10387 
.08 7819020 409516 6 1.0977378 12551 
.09 10250440 465725 -— 3 1.1138938 14372 
10 13147020 518852 —12 1.1314702 15237 
Formulae 


T = e*Ei(x) —+ 
S = xe Ei (x) 
Interpolation 
fn = (1 — )fo + fi + Ee Smo? + Ey5m* + Motyo! + Mitys' 
(See BAAS Table II, Auxiliary Table I, for definitions and values of E and M.) 
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EpitoriAL Note: It is easy to find the value of T(— x) for x positive 
and large, by use of the Laguerre approximate quadrature. We have 


T(-2) = eEi(-2) tte et — [ete tora. 


Now, approximately, 


(1) f (x 4. t)—e-* dt an Lax + x}? 
0 i=1 
where the x; are the zeros of the Laguerre polynomial and the a; are 
the corresponding Christoffel numbers. These have been tabulated by 
SALZER & ZUCKER.‘ 
The error in the formula (1) above is the value of 
_ {nid —* tint 
E, = (2n)! de" (x + t) — (n!) (x + t) 


for some i, 0 < t < «. This error is certainly less than 
E, = (n!)2/x29+1 
and since 
Enys/E, = (n + 1)*x 
the best ” to use is m = [x] which gives 
E,~ 2ne-™ 


which is about 10-* in the case x = 10. The accuracy obtained is therefore 
comparable with that of the above table for z = .1 and z = .09 and better 
when z > .08. 

This method will be as convenient in the evaluation of functions of the 
same general form as Ei(—x) even when the differential equation which 
they satisfy is complicated, or awkward to obtain. It cannot, however, be 
applied conveniently to the case of negative x. 

i% x 


The es of this note was sponsored in part by the Office of Air Research, 
AMC, USAF. 


 « R. Arrey, “The converging factor in asymptotic series and the calculation of 
Bessel, Laguerre and other functions,” Phil. M Mag., v. 24, 1937, p. 521-552. 

2A. C. Coutson & W. E. DuNcANsON, “Some new values for the exponential integral,” 
Camb. Phil. Soc., ey v. 33, 1942, . 154-761. 

*L. Fox & E. T. ‘Goopwin, “Some new methods for the numerical integration of 
ordinary differential equations,’ * Camb. Phil. Soc., Proc., v. 45, 1949, p. 373-388. 

4H. E. Satzer & RutH Zucker, “Table of the zeros and weight factors of the first 
fifteen Laguerre polynomials,” Amer. Math. Soc., Bull., v. 55, 1949, p. 1004-1012. 

5 G. SzEG6, Orthogonal Polynomials. Amer. Math. Soc., Coll. Publ., v. 23, 1939, p. 369. 


BIBLIOGRAPHY Z-XVI 


1. Anon., “Standards on electronic computers: definitions of terms, 
1950,” I.R.E., Proc., v. 39, March 1951, p. 271-277. 


2. Epmunp C. BERKELEY, “The uses of automatic computers in financial 
and accounting operations,’ Journal of Accountancy, v. 90, October 
1950, p. 306-311. 


The author feels that the needs of business in accounting are beginning 
to create a demand for the large electronic computers such as are now being 
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used in scientific and engineering problems. He explains, in this article, the 
fundamental nature of the large computer by comparing it in its essential 
characteristics with the ordinary desk calculator. 

By analyzing the basic requirements of typical accounting problems, 
he is able to forecast the changes in computer design necessary to meet 
these requirements. The tangible benefits to business, such as the lower 
cost of the information produced and quick access to completely assimilated 
up-to-the-minute information, are forcibly pointed out. Of course, there will 
be construction and training delays, and it would be wise for business to 
study some of the problems suggested in the paper so that the machines 
can go to work as effectively as possible after their construction. 

J. Blum 
NBSCL 


3. ENGINEERING RESEARCH ASSOCIATES, Inc., High Speed Computing 
Devices. McGraw-Hill, 1950, 451 pages. 


As stated in the preface, this book is primarily a discussion of the 
mechanical devices and electrical circuits which can be incorporated into 
computing machines. Included also are a brief treatment of the concepts 
of number theory which are pertinent to many electronic digital computing 
machines now in use or under development, a chapter on numerical analysis 
listing methods for solving algebraic and differential equations, and several 
chapters on computing systems. 

The book contains three parts: Part I, The Basic Elements of Machine 
Computation; Part II, Computing Systems; and Part III, Physical Com- 
ponents and Methods. It appears to have been directed mainly at readers 
with some knowledge of engineering who wish to learn how large-scale 
digital computing machines function, and it serves this purpose well. The 
treatment of analogue devices is quite cursory, however, and the modicum 
of numerical analysis included serves best to justify the definition of com- 
ponent, given in the introduction, as ‘‘any physical mechanism or mathe- 
matical method which is used as a tool in automatic computation.” 

Granting the definition of component as given in the book, the classi- 
fication of the material contained therein which is made in the introduction 
appeals to this reviewer. This classification of contents is the following: 

A. General. Chapter 1, Introduction, Chapter 2, Preliminary Considera- 
tions; B. Mathematical Components. Chapter 6, Arithmetic Systems, 
Chapter 7, Numerical Analysis; C. Physical Components and Methods. 
Chapter 3, Counters as Elementary Components, Chapter 4, Switches and 
Gates, Chapter 13, Arithmetic Elements, Chapter 14, Transfer Mediums, 
Chapter 15, Data Conversion Equipment, Chapter 16, Special Techniques 
and Equipment for Possible Use in Computing Systems, Chapter 17, Factors 
Affecting Choice of Equipment; D. Computing Systems. Chapter 5, A 
Functional Approach to Machine Design, Chapter 8, Desk Calculators, 
Chapter 9, Punched-card Computing Systems, Chapter 10, Large-scale 
Digital Computing Systems, Chapter 11, Analog Computing Systems, 
Chapter 12, The Form of a Digital Computer. 

As the chapter titles indicate, the book by the staff of the Engineering 
Research Associates touches but lightly on control and input-output ele- 
ments of high-speed computing devices. By way of compensation, however, 
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the arithmetic element of such devices is treated in detail and Chapter 14 
on transfer mediums, the largest chapter in the book, provides an excellent 
description of the physical phenomena basic to the high-speed storage 
elements in use in large-scale electronic computing machines. To those who 
are concerned most with electronic computing devices, Chapter 14 is prob- 
ably the most interesting chapter of the book, covering both storage elements 
of the dynamic type, like the acoustical delay line and the magnetic drum, 
and electrostatic storage elements depending, for their action, on the 
secondary-emission phenomenon. The book, prepared by competent staff 
members of a company active in the development of large-scale, high-speed 
computers, has real merit in more than a single chapter, however. In the 
opinion of this reviewer, it is a valuable addition to the literature on com- 
puting devices, which will be most useful to those interested in large-scale 
electronic digital computers. 
E. W. C. 


4. S. Gu, “A process for the step-by-step integration of differential 
equations in an automatic digital computing machine,’’ Camb. Phil. 
Soc., Proc., v. 47, 1950, p. 96-108. 


As the author observes, many of the methods of integrating differential 
equations which are effective for hand computation have serious disadvan- 
tages when the integration is to be performed on a high-speed automatic 
machine. Among the disadvantages may be listed: 


(1) Starting values may be obtained by a process different from that 
emphasized in the remainder of the integration. The extra orders required 
to do this on a machine are wasteful of memory space. 

(2) Preceding functional values may be required for the calculation at 
each point. The transferring of values and the “searching’’ which this 
necessitates is cumbersome and expensive of orders and time. 

(3) Changing the interval may be complicated. Here again, the number 
of orders required to effect the change may be excessive. 


These problems are particularly vexing when we are seeking to solve a 
large system of differential equations. The author, by an application of 
processes developed by Kutta, has worked out an extremely effective 
method, having none of the above disadvantages, by which such systems 
can be solved. It is straightforward, easily coded and economical of storage 
space. The paper includes a fairly thorough discussion of the accuracy 
obtainable in such a scheme. 

M. MONTALBANO 
NBSCL 


5. D. R. HartrREE, “‘Automatic calculating machines,” Math. Gazette, 
v. 34, 1950, p. 241-252, illustrs. 


This paper is the text of an address given before the Mathematical 
Association. Because this discussion was given before a mathematical group, 
the treatment of automatic calculators in general was limited. The author 
does touch upon some of the more important aspects of automatic machines 
with special reference to BABBAGE’s dream of an ‘‘Analytical Engine,” and 
the first realization of that dream, the Harvard Mark I Calculator. A short 
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discussion of the first machine using electronic circuits, the ENIAC, with 
illustrations of the main parts, completes the first part of the paper. 

The second part is concerned with a discussion of the EDSAC, which 
is a serial binary machine with a storage system using ultrasonic acoustic 
delay units. This machine uses one-address code instructions and serial 
storage. It has one storage unit for current instructions and another for 
registering the address from which the current instruction is taken. Nor- 
mally the contents of the latter are increased by unity for every instruction 
carried out. Numbers and instructions are received by the machine in coded 
decimal form by means of teleprinter tape, and output is on a teletypewriter. 
The numbers are operated on in the machine in the binary system. 

The author concludes with a note on programming and coding and 
gives an example of coding for the EDSAC. 

In this reviewer's opinion, the closing comment of the author is worthy 
of quoting. ‘““These machines have been called ‘electronic brains,’ which 
carries the suggestion that they can ‘think for themselves,’ which they 
cannot do: they can only carry out, quite blindly, the sequence of instruc- 
tions which has been thought out for them.” 


B. F. HAnpy, Jr. 
NBSCL 


6. MARSHALL KINCAID, JOHN M. ALDEN, & RoBerT B. Hanna, “Static 


magnetic memory for low-cost computers,”’ Electronics, v. 24, 1950, 
p. 108-111. 


This paper reviews the basic theory and gives several applications for 
the static magnetic storage units using the special magnetic material Delta- 
max. The authors feel that these magnetic storage units have application 
in medium-speed, low-cost computers which could handle many problems 
encountered in daily operation of business and industry. The applications 
cited in the paper for use of these devices in computers are as follows: 


1. Storage of information. The handling rate varies from zero to 30,000 
pulses per second, and the informaticn is not lost in the event of power 
failure. 

2. Transformation of information. The devices transform information 
slow-speed to high-speed pulse systems and vice versa; they also change 
parallel-type information to serial-type information and vice versa. 

3. A high-speed counter which can count up to 25,000 pulses per second 
using only simple auxiliary circuitry. 


The paper concludes by mentioning some of the production techniques 
and difficulties encountered in the commercial production of these devices. 


M. M. ANDREW 
NBSMDL 


7. B. L. Moore, “Pentode counting or control ring,” Review Sci. Instr., 
v. 21, 1950, p. 337-338. 


A counting or control ring circuit using 6AS6 pentodes in a very stable 
manner is described. The input signal is said to be essentially independent 
of amplitude and wave shape, and the values of circuit components and 
supply voltage are not critical. Control voltages may be taken from one or 
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more stable states of the counter, the maximum number being limited 
essentially by the power required to drive the grid. Combinations of counters 
may be used to obtain a ring of many stable states. 


M. M. ANDREW 
NBSMDL 


8. Ropert F. SHaw, “Arithmetic operations in a binary computer,” 
Review Sci. Instr., v. 21, 1950, p. 687-693. 


The author gives a summary of the methods used to facilitate the basic 
arithmetic computations in contemporary binary computers. The use of 
complements is explained in defining negative numbers for machine use; 
corrections for multiplication error and round off procedures are discussed. 
With numerical examples, simply written but thorough, this paper is a good 
introduction to the logical mechanism of the binary computer. 


KARL GOLDBERG 
NBSCL 


9. M. V. WILKEs, “‘Automatic Computing,” Nature, v. 166, December 2, 
1950, p. 942-944. 


Describes a Summer school course in coding for an automatic digital 
computing machine held at Cambridge in the University Mathematical 
Laboratory from September 12th through 21st. The course was concerned 
mainly with the methods used in connection with the EDSAC. 


NEws 


Conference on automatic computing machinery and applications.—The conference 
was jointly sponsored by the Association for Computing Machinery, the Industrial Mathe- 
matics Society, and the Advisory Committee for Wayne University Computation Labora- 
tory. There were general sessions on Tuesday morning and evening, March 27th, and two 
parallel groups of sectional meetings beginning Tuesday afternoon and continuing through 
Wednesday. 

As the conference was held in the industrial center of Detroit, there was considerable 
interest in the applicability of computing machines to industrial problems. An exhibit of 
computer equipment was shown Monday evening through Wednesday afternoon. The 
program for the meeting was as follows: 





General Session, Tuesday, March 27 
Welcoming address 


“Some computing problems in the automo- 
tive industry” 

‘Automatic calculations and their applica- 
tions”’ 

Sectional Meeting, Tuesday, March 27, Sec- 
tion A 

“Digital computer research at M.I.T.” 

“The ENIAC—a five year operating survey” 


“Problems solved on B.T.L. Model VI Com- 
puter”’ 

“Capabilities and limitations of existing 
electronic equipment” 

Sectional Meeting, Tuesday, March 27, Sec- 
tion B 


Franz L. Att, President, ACM, Chairman 

Victor A. Rapport, Dean, Wayne Uni- 
versity 

Paut T. Nims, Chrysler Corporation, Vice- 
president, Industrial Mathematics Society 

Howarp H. AIKken, Director, Computation 
Laboratory, Harvard University 

Avex L. Haynes, Ford Motor Co., Chair- 
man 

Rosert R. Everett, M.I1.T. 

W. BarKLEy Fritz, BRL, Aberdeen Proving 
Ground 

Ernest G. ANDREWS, Bell Telephone Lab- 
oratories 

Ipa R. Ruopes, NBS 


C. C. Hurp, IBM, Chairman 





“‘Floating decimal calculations on the card 
programmed electric calculator” 

‘The treatment of systems of differential 
equations on the SEAC” 

“Solution of the non-linear supersonic flow 
equations on the SEAC” 

“Preparation of a theodolite problem for 
large-scale machines” 

Evening Session Dinner, Tuesday, March 27 


“Electronics for business—luxury or neces- 
sity?” 

Sectional Meeting, Wednesday Morning, 
March 28, Section A 


“The programming of some matrix com- 
putations for the Mark II Aiken Relay 
Calculator” 

“The operation of the Fairchild specialized 
digital computers for certain applications 
in matrix algebra” 

‘Matrix algebra programs for the UNIVAC” 


“A UNIVAC program for inventory require- 
ments” 

Sectional Meeting, Wednesday Morning, 
March 28, Section B 

“Automatic checking features of the Ray- 
theon Digital Computer” 

“The external memory of the Raytheon 
Digital Computer” , 

“The digital reader” 

“Preliminary considerations on a magnetic 
drum controlled computer”’ 

“Techniques in the design of digital com- 
puters” 

Sectional Meeting, Wednesday Afternoon, 
March 28, Section A 

“‘On the accuracy of Runge-Kutta’s method” 


“Optimum trajectories” 
“Reversing digit number system” 


‘Analysis of digital computers in control 
systems” 

Sectional Meeting, Wednesday Afternoon, 
March 28, Section B 

**The impact of business on computing ma- 
chinery design” 

‘*Plans for tabulating parts of the 1950 popu- 
lation census on an electronic computer” 

“Engineering applications on the 
MADDIDA” 

“Mathematical techniques in the uses of 
MADDIDA” 
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Ricwarp H. Starx, Los Alamos Scientific 
Lab. 
Joseru H. Levin, NBS 


ErHet C. Marpen, NBS, 

J. Conrap Crown, NOL 

BERNARD DimsDALE, BRL, Aberdeen Prov- 
ing Ground 

Joun R. RicHarps, Wayne University, 
Toastmaster 

Joun S. CoLeman, President, Burroughs 
Adding Machine Co. 

ROBERT SCHILLING, General Motors Corpo- 
ration, President, Industrial Mathematics 
Society 

ALLEN V. HERsHEy, Naval Proving Ground 


J. J. Stone, Fairchild Engineering and Air- 
plane Corporation 


H. F. MitcHe tt, Jr., Eckert-Mauchly Com- 
puter Corporation (Subsidiary, Remington 
Rand Corp.) 

G. M. Hopper, Eckert-Mauchly Computer 
Corp. Subsidiary, Remington Rand Corp.) 

ALBERT C. HALL, Bendix Aviation Corpo- 
ration, Chairman 

Louts Fern, Raytheon Manufacturing Com- 
pany 

KENNETH M. REBLER, Raytheon Manufac- 
turing Company 

GritBert W. Kino, Arthur D. Little, Inc. 

THEODORE SHAPIN, Jr., E. F. Moore, Uni- 
versity of Illinois 

R. E. Spracue, Computer Research Cor- 
poration 

C. C. Bramsie, Naval Proving Ground, 
Chairman 

Max Lorkin, BRL, Aberdeen Proving 
Ground 

ARNOLD S. MENGEL, The Ranp Corporation 

GerorGE W. Patterson, Burroughs Adding 
Machine Company 

W. K. Lrnvitt, J. M. Sauzer, M.I.T. 


Epmunp C. BERKELEY, Secretary, Associa- 
tion for Computing Machinery 

E. F. Coo.ey, Prudential Insurance Com- 
pany of America 

James L. McPaerson, Bureau of the Census 


GLENN E. HacGen, Northrop Aircraft, Inc. 


Myron J. MENDELSON, Northrop Aircraft, 
Inc. 
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Institute on Computing Methods and Machines.—On January 26th through 30th 
inclusive a duplicate meeting was held at the University of California at Berkeley (Jan. 
26-27) and at Los Angeles (Jan. 29-30). The purpose of the meeting was to provide a 
general picture of computer application to problem solutions and a specific picture of the 
computer facilities available in the University of California. The meeting was co-sponsored 
by various departments of the universities and by the NBSINA. The following program 
was presented: 


Program Berkeley Los Angeles 
Morning Session 
(Jan. 26 at Berkeley and Jan. 29 at Los 
Angeles) 
Chairman H. A. Scuape, U.C., L. Caerry, NAMTC, 


Welcome and Introduction 


Berkeley 

M. P. O’Brien, 
Chairman, Dept. of 
Eng., U. C., Berkeley 


Point Mugu, Calif. 

L. M. K. Boe.tTer, 
Chairman, Dept. of 
Eng., UCLA 





“Current status of digital computer de- E. U. Connon, 
velopments” Director, NBS 
“The facilities and program of the Com- Paut L. Morton, 
puter Laboratory at Berkeley” Director, Computer 
Lab., U.C., Berkeley 


“The California Digital Computer— Davip R. Brown, 


CALDIC” U.C., Berkeley 
Afternoon Session 
Chairman Jerzy NEYMAN, Louts A. Prres, UCLA 


Director, Statistical 
Lab. U. C., Berkeley 
“Role of the mathematics department in Derrrick H. LEHMER, 
a university computing center” U.C., Berkeley 
“Computations relating to some solar LELAND E. CUNNINGHAM, 


system problems” U.C., Berkeley 
“Facilities of the Institute for Numerical 
Analysis: 
a) problem formulating service A. S. Cann, 
NBSINA 
b) punched card equipment E. C. YowE i, 
NBSINA 
c) electronic equipment—SWAC” Harry D. Huskey, 
NBSINA 
Inspection of Computer Facilities 
Morning Session 
(Jan. 27 at Berkeley and Jan. 30 at Los 
Angeles) 
Chairman ALFREDO Banos, Jr., WitiiaM A. MERSMAN, 
UCLA Ames Aeronautical 
Lab. 
“Recent developments in analog com- W. W. Soroka, 
puter techniques” U.C., Berkeley 
“The network thermal analyzers” F. W. Scuortrt, 
UCLA 
“Electronic analog computations” T. A. RoGERs, 
UCLA 
“Mechanical differential analyzer-prob- DoNnaLp LEBELL, 
lem solution” UCLA 


Inspection of Computer Facilities 
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The Institute of Radio Engineers.—At the National Convention of the Institute held 
in New York on March 19th through March 22nd, two sessions were devoted to computers. 
The programs were as follows: 


Computers I E. R. Prore, Office of Naval Research, 
Washington, D. C., Chairman 
“The Raytheon selection matrix for com- KENNETH M. REBLER, Raytheon Manufac- 


puter and switching applications” turing Company, Waltham, Mass. 
“Saturable reactors as substitutes forelectron James G. Mies, Engineering Research 
tubes in high-speed digital computers” Associates, Inc., St. Paul, Minn. 
‘Ferromagnetic cores for three-dimensional W1L1AM N. Paptian, Servomechanisms Lab- 
digital storage arrays” oratory, MIT, Cambridge, Mass. 
“‘Dependable small-scale digital computer” J. J. ConNnoLLy, Teleregister Corporation, 
New York, N. Y. 
‘*An asynchronous control for a digitalcom- D. H. Griptey, Naval Research Labora- 
puter” tories, Washington, D. C. 
Computers II J. W. Forrester, MIT, Chairman 
‘*A sampling analogue computer” JoHN BRooMALL and LEon RIEBMAN, Moore 
School of Electrical Engineering, Univ. of 
Penn. 
“A time division multiplier for a general- R. V. Baum and C. D. MorriLt, Goodyear 
purpose electronic differential analyzer” Aircraft Corp., Akron, Ohio 
‘*A high-speed product integrator” ALAN B. MacNEE, University of Michigan, 
Ann Arbor, Michigan 
“Plug-in units for digital computation” G. Gitnsk1 and S. Lazecx1, Computing 
Devices of Canada, Ltd., Ottawa, Ontario, 
Canada 
“A five-digit parallel coder tube”’ J. V. HARRINGTON and K. N. WULFSBERG, 


Air Force, Cambridge Research Labora- 
tories, Cambridge, Mass., and G. R. 
SPENCER, Philco Tube Laboratory, Lans- 
dale, Pa. 


National Physical Laboratory.—In November 1950 at the National Physical Labora- 
tory a pilot model of the ACE (Automatic Computing Engine) was demonstrated to 
members of the British Press. The ACE is an electronic digital computer which is being 
constructed under the aegis of the Department of Scientific and Industrial Research by 
the staff of Electronics and Mathematics Divisions of the National Physical Laboratory 
at Teddington, Middlesex, in association with the English Electric Company of Stafford. 

Owing to the size and complexity of the ACE itself, it was thought desirable to build 
first a smaller version known as the ACE pilot model. This is a serial machine operating 
on numbers in the binary scale, with a pulse.repetition rate of one megacycle per second. 
It uses about 800 tubes in 40 plug-in units mounted on a rack about 12’ X 16’ and has 
separate control panel and power supply units. The internal memory is of the mercury 
delay line type and consists of eight long tanks, each with a delay of 1024ysec. and 10 
short tanks, each with a delay of 32 usec. Storage capacity is 264 words, each word being 
the binary equivalent of a 10-decimal-digit number or an instruction in a three-address 
code. The time necessary to multiply two 10-decimal-digit numbers is .002 sec. For the 
input and output sections of the machine, modified Hollerith punched card equipment is 
used. The demonstration included the extraction of 3rd and 4th roots of 7-decimal-digit 
numbers, and a program for testing whether any given number less than a million is prime. 
Numbers, suggested by the audience, were fed into the machine in decimal form by hand 
switching. If the number was not prime, the lowest factor was indicated in decimal form 
by means of lights; if it was prime, after three seconds the machine sounded a buzzer. 
The machine was also used to solve a problem of practical importance—i.e., the tracing 
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of rays through a compound lens. The machine computed in a few minutes, the paths of 
20 rays through a compound lens with six refracting surfaces. Such a computation would 
have taken a human computer two 7-hour days with a desk calculator. 


Reeves Instrument Corporation.—Project Cyclone Symposium I on REAC Techniques 
was held in New York City, March 15-16, 1951 under the sponsorship of the Reeves 
Instrument Corporation with the approval of the U. S. Navy Special Devices Center. 

It was the purpose of Project Cyclone Symposium I on REAC Techniques to provide 
a means for discussing techniques, applications and engineering difficulties which have 
arisen among users of REAC equipment, and to stimulate new techniques in the field of 
analogue computation. 

The program consisted of five sessions, under the chairmanship of Stanley Fifer, head 
of Project Cyclone. 


Morning Session, Thursday, March 15: 


Welcoming addresses by David T. Bonner, President, Reeves Instrument Corporation, 
and Paul Staderman, Head, Synthetic.Warfare Section, Special Devices Center. 


1. ‘Optimum trajectories,”” Arnold S. Mengel, The RAND Corporation. 

2. “Solution of some partial differential equations on the REAC,” G. W. Evans II, 
Argonne National Laboratory. 

3. “Analytical and mechanical methods of solution of differential equations with dis- 
continuous forcing functions,”’ William A. Mersman, National Advisory Committee 
for Aeronautics, Ames Aeronautical Laboratory. 

4. “The solution of polynomial equations on the REAC,” Louis Bauer and Stanley Fifer, 
Project Cyclone. 

5. “Computation of multidimensional integrals on the REAC,” Leslie C. Merrill, Argonne 
National Laboratory. 


Afternoon Session, Thursday, March 15: 


6. ““RAND REAC modifications,” Arnold S. Mengel and Wesley S. Melahn, The RAND 
Corporation. 

7. “‘Some modifications for maximum utilization of the REAC,” John L. Burnside, North 
American Aviation, Inc. 

8. ‘The use of the REAC by Minneapolis-Honeywell in the analysis and synthesis of 
automatic flight control,’’ Remus N. Bretoi and David L. Markusen, Minneapolis- 
Honeywell Regulator Company. 

9. “‘Slaving of a rotational mount to an analog computer,’ C. F. Cook, National Bureau 
of Standards. 

10. ‘Preliminary investigation of suitability of REAC for experimental curve fitting,” 
Donald S. Teague, Jr. and R. D. Gilpin, U. S. Naval Air Missile Test Center, Point 
Mugu. 

11. “‘The generation of an N dimensional normal distribution by means of analog equip- 
ment,” Harry H. Goode, William A. Wheatley and George G. den Broeder, Willow Run 
Research Center, University of Michigan. 


Morning Session, Friday, March 16: 


12. “The generation and measurement of ultra low frequency random noise,” Robert R. 
Bennett and A. S. Fulton, Hughes Aircraft Company. 

13. “Simulation of noise,’”’ Harold L. Ehlers and Erwin Vogel, Fairchild Engine and Air- 
plane Corporation. 

14. “Frequency analysis by electronic analog methods,” William A. McCool, Naval 
Research Laboratory. 

15. “Application of REAC equipment to the solution of problems involving integral 
operators, I,” Paul Brock and Seymour Sherman, Project Cyclone. 
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Noon Session, Friday, March 16: 
Numerical checks; introduction by Paul Brock, Project Cyclone. 
Afternoon Session, Friday, March 16: 


16. ‘‘A REAC solution of a linear proton accelerator design problem,’”’ A. H. Miller, pre- 
sented by Nathanial B. Nichols, University of Minnesota. 

17. “‘Representation of functions of several variables on REAC equipment,” Hans Meissinger, 
Project Cyclone. 

18. “‘A method for solving problems on the REAC by the use of transfer functions without 
passive networks,”’ Cyrus Beck, Naval Air Experimental Station, Philadelphia. 

19. ““REAC techniques,” W. Frank Richmond, Jr., The Glenn L. Martin Company. 

20. “‘The generation of straight line transfer relationships,” Robert R. Bennett, Hughes 
Aircraft Company. 


UNIVAC Acceptance Tests.—On March 23rd the first UNIVAC successfully passed 
its acceptance tests in the presence of representatives of the Bureau of the Census and the 
Bureau of Standards. A discussion of these tests will follow a brief description of this new 
computer. 

The UNIVAC system consists of a computer with an elaborate Supervisory Control 
Console, one to ten Uniservos, one or more Unitypers, and one or more Uniprinters. There 
is also available a Punch-Card-To-Tape Converter, and a device to perform the reverse 
conversion is under construction. The computer is a decimal machine containing 1000 
words in its 100 mercury delay lines. Each word consists of 12 alphanumeric digits and has 
room for two one-address commands. The pulse repetition rate of the machine is 2.25 
megacycles allowing about 2000 additions or over 450 multiplications to be performed in 
one second. Each half of the twin Arithmetic Unit checks the operations of the other half 
and indicates on the Supervisory Console any discrepancy that may occur. 

Each Uniservo uses an 8-channel metal tape with a magnetic plating. The tape is 
initially examined for flaws, which are marked and which the machine is constructed to 
skip over. These tapes move either forwards or backwards at ten feet per second, so that 
over 800 words can be delivered to the Memory or recorded from it every second. A reel 
of tape is over 1500 feet in length, has a capacity of 120,000 words, and may be reused 
some 700 times. It is completely erasable, moves in either direction, and can be rewound, 
all at the behest of the coder’s commands. 

The Unityper’s keyboard possesses a number of keys in addition to the conventional 
alphanumeric characters. One set of these keys serves to prepare commands leading to the 
desired format in which the final results are to appear; another set calls into play the group 
of commands which appears on three tracks of the rudimentary memory built into the 
Unitypers. These commands can be set to facilitate, shorten, and correct the typing of the 
input information. 

The Uniprinter converts the binary-coded information on the tapes into alphanumeric 
characters at the rate of 11 per second. At the same time it is able to carry out printing 
instructions stored on the tape such as: carriage return, upper case, ignore (i.e., do not 
print), etc. It should be added that a fast line printer is under construction which will list 
information at about twice the speed of the standard punch-card tabulator. 

The Card-to-Tape Converter handles six 80-column cards per second, storing the 
contents of 10,000 cards on a reel of tape and stopping automatically after a preset number 
of cards have been converted. 

The four acceptance tests which the UNIVAC system underwent successfully con- 
sisted of: 


1. A general test with emphasis on the computer’s performance. It was subdivided 
into three parts. The first, which was repeated 404 times per test, investigated the 
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reliability of the execution of every one of the computer’s commands. The second 
part required the computer to sort a number of items obtained as a result of the 
operations in the first part, as well as to solve a rather complicated partial differential 
equation. In the last part, the tapes of the four Uniservos were run through a difficult 
routine of a series of movements, delivering and recording information. Two such 
tests, requiring twenty minutes for a flawless performance, were designated as a 
unit. Out of the successive nineteen units which the UNIVAC performed during 
the test period, sixteen were carried out perfectly; in the three units where stops 
occurred due to the action of the error-detection circuits, the operator was able to 
rectify them by the mere manipulation of the console buttons. 

2. A Uniprinter test requiring a sentence and a numerical table to be printed out in 
proper form 200 times. This test lasted eight hours, during which five errors occurred, 
whereas the conditions of the test allowed eight. 

3. A Card-to-Tape-Converter test, which was passed satisfactorily. Although required 
to process a deck of 10,000 very badly punched cards, the Converter made only one 

mistake while storing this information repeatedly on ten tape reels. : 

. A general test, with emphasis on Uniservo performance. Various aspects of tape 

movements were under scrutiny such as any tendency of the tapes to move out of 
position during prolonged inactivity, the efficiency of flaw-detection in the magnetic 
coating of the tapes, the resistance of this coating to a continuous rerunning of the 
tape under the reading heads for 700 times, the effect that various juxtapositions 
of tape commands might have on their performance, and several other points of 
possible weakness. This test lasted ten hours during which six errors were detected 
by the machine. 


de 


The machine will remain on the premises of the Eckert-Mauchly Computer Corporation 
for about a year, performing computations for the Bureau of the Census and operating 
24 hours a day. 


OTHER AIDS TO COMPUTATION 
BIBLIOGRAPHY Z-XVI 


10. B. P. BoGert, “‘A network to represent the inner ear,’’ Bell Lab. Record, 
v. 28, 1950, p. 481-485. 


An electrical network consisting of 175 sections is used to simulate the 
cochlea which is believed to be the frequency sensing portions of the ear. 
The results of measurements on this simulator are compared for a number 
of frequencies with the results computed from the theory developed in a 
paper by PETERSON & Bocert.' The simplified model upon which a hydro- 
dynamical theory is based consists of two parallel ducts separated by a 
movable membrane of varying width. Measurements on the electrical net- 
work are compared with the numerical solution of the equations of wave 
propagation in this model, in regard to the amplitude of the response of the 
membrane at different points along its length to different frequencies. The 
response patterns vary considerably with the frequency but good agreement 
was obtained between the network, the theoretical model and certain 
experimental results for the actual cochlea. 


F. J. M. 


1L. C. Peterson & B. P. Bocert, “A dynamical theory of the cochlea,” Acoust, 
Soc. Amer., Jn., v. 22, 1950, p. 369-381. 
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11. J. W. CrarK & R. E. NEusBeEr, “A dynamic electron trajectory tracer,” 
I.R.E., Proc., v. 38, 1950, p. 521-524. 


In the gravitational analog of a vacuum tube, a surface reproduces the 
electric potential in the tube and electrons are simulated by metal balls. 
By vibrating this surface, the authors reproduce the effect of the application 
of radio frequency voltages to various electrodes. Relatively low mechanical 
frequencies may be used. The results are recorded by a stroboscopic photo- 
graphic process. A number of such photographs are given as illustrations. 

F. J. M. 


12. H. L. Curtis, ‘Determination of curvature by an osculometer,’”’ NBS, 
Jn. Research, v. 44, 1950, p. 131-134. 


The radius of curvature of a curve is determined graphically by means of 
a series of pairs of arcs of known curvature drawn on transparent paper. The 
user tries to arrange a pair of arcs symmetrically on both sides of the given 
curve. “The curvature at a point can usually be determined to 10-* recip- 
rocal centimeters, provided the curvature is approximately constant for a 
distance of 2 cm. . . .”” A discussion is given for the graphical determina- 
tion of accelerations by this method. 

F. J. M. 


13. T. S. Gray & H. B. Frey, ‘Acorn diode has logarithmic range of 10°,” 
Rev. Sci. Inst., v. 22, 1951, p. 117-118. 


The diode type 9004 with a d.c. filament voltage supply of 4 volts has a 
logarithmic current characteristic between 10-* and 10- amperes for the 
plate voltage range of 0 to —1.75 volts. 

F. J. M. 


14. H. Iams, “A method of simulating propagation problems,” I.R.E., 
Proc., v. 38, 1950, p. 543-545. 


Propagation problems can be effectively simulated by propagating 
centimeter wavelength radio waves in large sheets of low loss dielectric 
materials. The conditions of the problem are simulated by imbedding 
objects in the sheet or utilizing the effect of varying the thickness of the 
sheet to simulate a change in refractive index. The propagation is recorded 
by means of a “phase front plotter’’ which records the lines of equal phase. 
A coupling set up to simulate waves propagated from a large distance is 
described. 


F. J. M. 


15. G. A. Korn, “Stabilization of simultaneous linear equation solvers,” 
I.R.E., Proc., v. 37, 1949, p. 1000-1002. 


The author gives an experimental test involving only one equation on 
an amplifier such that stability in this test will permit the amplifier to be 
used in a linear equation solver of the Goldberg-Brown type for any number 
of unknowns. A discussion of this paper, by L. A. ZADEH and the author, 
will be found in I.R.E., Proc., v. 38, 1950, p. 514. 


F. J. M. 
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16. G. LrEBMANN, ‘Solution of partial differential equations with a resist- 
ance network analogue,” Brit. Jn. Appl. Phys., v. 1, 1950, p. 92-103. 


The author describes a resistance network analogue which can be used 
for the solution of Laplace’s equation. The author points out that such a 
network may be considered as a computing machine for carrying out the 
finite difference method. The expected accuracy of the resistance network 
solutions is influenced by the finite mesh size and by resistors which deviate 
from the theoretical value. The error caused by finite mesh size can be 
greatly reduced by an extrapolating method due to Richardson. This 
method is briefly discussed in the paper. An examination of the errors due 
to deviations in resistance values follows but in order to evaluate this error 
the author must assume random distribution of deviations from the theoreti- 
cal value. With this assumption he finds that the errors are exceedingly 
small. However, since the distribution of resistance deviations cannot 
usually be predicted definitely, the method of the author does not permit 
the determination of the maximum error incurred with a network. 

In describing the design of a network he presents a valuable method of 
overcoming the difficulty of transition from coarser to finer lumping (his 
Figure 11). The actual description of a network with which the author 
works shows that he has assembled more than 2600 resistors. The descrip- 
tion is all too brief and gives the impression that the author has designed 
a network with fixed resistance values instead of using as frequently done 
resistors which could be set to desired values. Using fixed resistors would 
seriously limit the applicability of his apparatus. He uses alternating current 
as power supply and an oscilloscope for balance indicating. It would appear 
that direct current would be experimentally easier. The bibliography over- 
looks an important British publication which obviates the necessity in the 
case of a study of cylinders to replace a solid cylinder by a hollow one.' 
Moreover, he disregards the considerable literature on such networks used 
in this country. The Heat and Mass Flow Analyzer Laboratory at Columbia 
University has issued over the last 10 years some 50 papers all dealing 
with work on a resistance capacitance network which of course is a more 
inclusive system than the one dealt with in the paper under review. A 
reasonably complete bibliography of this laboratory up to 1947 is contained 
in reference.? 

V. PascHKIS 
Columbia University 
New York, N. Y. 


1R. Jackson, R. J. SarGant, J. B. Wacstarr, N. R. Eyres, D. R. Hartree & J. 
INGHAM, “Variable heat flow in steel,”’ Iron & Steel Inst., Jn., v. 150, 1944, No. 2, p. 211- 
268 


2V. Pascukis, “The heat and mass flow analyzer—tool for heat research,” Metal 
Progress, v. 52, 1947, p. 813-818. 


17. W. S. McEwan & S. SKoLNIK, “An analog computer for flame gas 
composition,” Rev. Sci. Inst., v. 22, 1951, p. 125-132. 


A problem in chemical equilibrium reduces to the problem of solving a 
system of equations, some of which are linear in the unknowns, others are 
ratios which may involve square roots of the unknowns. The device de- 
scribed realizes the unknowns as resistances and the equations are set up 
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as Wheatstone bridge circuits. While ganged potentiometers are used in 
certain instances, a switching system permits one to realize various different 
equations with the same resistance as unknowns in most cases. The device 
is manipulated until all the bridge circuits are in equilibrium. The manipu- 
lation is guided by certain physical information. 

F. J. M. 


18. A. Many, “An improved electrical network for determining the eigen- 
values and eigenvectors of a real symmetric matrix,” Rev. Sci. Inst., 
v. 21, 1950, p. 972-974. 


An electrical network for the determination of the characteristic roots 
and vectors of a symmetric matrix of order ten, consisting of ten LC-circuits 
coupled to each other by equal condensers, is described. (The characteristic 
roots correspond to the resonant frequencies of the LC-circuits.) A similar 
but smaller and less accurate machine was described by Many & MErBoom.! 
In the present paper the effects of unequal self-inductance of the coils of 
the tank circuits is compensated for. The stray capacitance of the wiring 
of the network is also compensated for by trimming the coupling condensers. 
Two important factors affecting the accuracy of the eigenvalues are: 
(1) Accuracy with which the matrix is realized, and (2) Losses in the coils. 
The above discussion compensates for (1). Using some results of the earlier 
paper, it is pointed out that the accuracy for matrices with reasonably 
well separated roots is mainly determined by (1) while for matrices with 
close eigenvalues, (2) predominates. In the first case, the accuracy obtained 
was about 0.01% of the difference of the greatest and smallest roots. 


K. S. MILLER 
New York, N. Y. 


1A. Many & S. MErBoom, “Electrical network for determining the eigenvalues and 
eigenvectors of a real symmetric matrix,” Rev. Sci. Inst., v. 18, 1947, p. 831-836. 


19. R. W. MARSHALL, “‘An integraph for semicurvilinear coordinate paper,” 
Bell Lab. Record, v. 28, 1950, p. 50-52, 65. 


When a function y(#) is obtained from the output of a customary re- 
cording instrument, it is plotted on curvilinear paper, i.e., the output is 
measured along an arc of a circle of fixed radius. The device described in 
this article has an input arrangement which changes this displacement to a 
linear displacement. This linear displacement is fed into a Corap1 integraph 
which plots a curve corresponding to the integral of the original function. 

F. J. M. 


20. G. E. REYNoLDs, A New Method for Square Root Evaluation, Particularly 
to Eight Significant Figures with Monromatic Calculator. Publication 
E4077 Unclassified. 22 p. Air Force Cambridge Research Laboratories, 
Cambridge, Mass., 1951. 21.4 X 27.6 cm. 


The booklet’s introduction states that it stems from ideas presented 
by R. C. Spencer in unpublished memoranda of 1948, and also that with 
some changes the method is adaptable to the Fridén and Marchant calcu- 
lators, as outlined in an appendix. 
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By use of a compact optimum-interval table of certain squares of selected 
numbers as well as of 4-times the numbers (in an adjacent column), a root 
with error not exceeding 5 in 9th place is obtained by comparatively simple 
calculator work. This results from combining two Newton-type iterations 
to a square root; i.e., the mean of quotient and divisor after dividing a num- 
ber whose root is desired N by an approximate root A obtains a root sub- 
stantially correct to twice the number of digits that are correct in the first 
approximation of the root. This is the first time so far as the reviewer 
knows that two successive Newton iterations have been combined for calcu- 
lator usage as a continuous process. Two divisions are required, as is to be 
expected. 

The optimum-interval table used has 138 lines of which the following 
are typical: 


APPROX. DIV. 
11.8336 13.76 
12.2500 - 14.00 
12.6736 14.24 


The left-hand number is the square of a first approximation of N?, or A?, 
and the right-hand number is 4A, which are applied as follows: 


(N + A*)/4A + 0.25N/((N + A?)/4A) = Ni+e 


The selection of A? and 4A values for the table is so made that ¢ does 
not exceed 5 in 9th place of N}. 

The approximation amount A? may be either the tabular value above or 
below N when the latter is pointed off in the manner of the tabular values. 

The method of obtaining the tabular intervals so as to minimize the 
size of the table is given, as are curves from which the intervals may be 
determined for preparation of tables for obtaining 7-, 8-, or 9-figure roots, 
though only the table for 8-figure roots is given. A feature of this part of 
the explanation is that a complicated expression involving the interval 
between tabular values and the approximate root for various table ranges 
reduces to an exceedingly simple form. The final formula used for calculator 
work, as above, is notable because the quotient of the first division becomes 
the divisor of the second. In the case of the Monromatic Calculator, a single 
entry of N suffices for both terms, it being held in the storage register for 
multiplication by 0.25 to obtain the numerator of the second term. 

The appendix states that the method is applicable to late-model Fridéns 
without change of formula, but a change with an additional step is recom- 
mended for the Marchant. The reviewer submits that a change is not 
required in this instance because even though the first quotient remains 
in the register so the second quotient may be accumulated, the keyboard 
entry of N for starting the second term may still be multiplied by 0.25 to 
obtain 0.25N in the product register, provided N is then cleared from the 
keyboard and a negative entry of 0.25 is made in the multiplier keyboard 
to restore the original quotient. This is readily done on any Marchant 
having automatic multiplication. 

The proposed machine method requires that certain exceptions have to 
be made to the general rule that the first significant figure of all entries is 
to be made at the extreme left of keyboard; no decimal markers apparently 
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are used. The reviewer suggests that unless extreme-left entries can be 
made without exception, it is better practice and requires less attention by 
the operator if there is a suitable decimal setup, and all entries are made 
with reference to it. As the table shows its amounts pointed off in the range 
1-100, any N may similarly be pointed off without loss of generality; i.e., 
82547.355 would be considered as 8.2547355. A suitable decimal setup then 
would be: counting register, at left of 8; product register, at left of 16; and 
keyboard, at left of 8; set No. 9 tab key. This setting accommodates at 
least 9-figure N’s and roots to 9 digits as well as intermediate amounts 
during the solution; an 11th dial in the counting register is not required. 
T. W. Smmpson 
66 Alvarado Road 
Berkeley, Calif. 


21. R. Roscog, ‘Mechanical models for the representation of visco-elastic 
properties,” Brit. Jn. Appl. Phys., v. 1, 1950, p. 171-173. 


The author establishes the existence of two canonical forms for the 
networks of springs and dashpots used to represent visco-elastic properties. 
One form is parallel, the other is essentially a series arrangement of parallel 
pairs of springs and dashboards with a simple series addition. The author’s 
discussion is based on the analogy with resistor-capacity circuits for which 
the corresponding result is known. (Cf. Guillemin, ‘Communication Net- 
works,”’ v. 2, Chap. V, p. 184-221, John Wiley and Son, New York, 1947.) 

F. J. M. 


22. L. H. Witson & A. J. MILEs, ‘Application of the membrane analogy 
to the solution of heat induction problems,” Jn. Appl. Phys., v. 21, 
1950, p. 532-535. 


The authors describe the application of the membrane analogy to the 
solution of heat conduction problems. Mathematical justification is obvious 
from the identity of the equations and is dealt with only briefly. In order 
to carry out the technique, a membrane is created by using a soap film. 
The deviation of distance of various parts of the membrane from a reference 
point is proportional to the deviation of temperature in an equally shaped 
plane from the temperature of the corresponding reference point. The 
membrane is stretched across appropriate boundaries; for example, if heat 
flow through a hollow cylinder is to be studied, the two surfaces of which 
are at different temperatures, an annular membrane must be stretched 
between two cylindrical bodies serving as boundaries and placed at different 
levels according to the different temperatures of the outside and inside 
surface of the cylinder. The authors show two examples of comparison of 
determinations carried out by this method against theoretical curves. In 
one case the check is very good and seems to be better than 1%. The other 
example shows deviations which are, in part, as large as 10%. 

Whereas the technique undoubtedly has some merits, the following 
serious limitations are not discussed and not even mentioned in the paper. 
First, in many thermal problems, the surface temperatures are not given; 
instead the temperature of the ambient is known and a boundary conduct- 
ance applies, which may or may not be introduced as a function of the surface 
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temperature. Application of such boundary conductances seems impossible 
with the present technique. Secondly, in many steady-state problems two 
or more materials of different conductivity enter the picture. Application 
of the membrane technique to such cases does not seem possible. The 
authors mention the desirability of investigating the membrane analogy in 
the realm of transient-heat conduction problems. There is no indication 
either on a mathematical or experimental ground that such approach is 
possible. 
V. PAscHKIS 

Columbia University 

New York 27, N. Y. 

NOTES 


127. AN ITERATIVE Process.—lIn a recent paper HARTREE' calls atten- 
tion to the need for the development of a general method of obtaining 
iteration formulas of second and higher orders for the solution of any 
“algebraic” (i.e., not differential) equation. If y = Y is the true solution 
of F(y) = 0, and y, = Y +, is an approximation to Y, the formula 
¥n41 = G(ya) is an iterative process. If a4: = O(n,"), Hartree calls the 
process r-th order. The number of correct figures in y,4; is approximately 
r times the number in y,, so it is obvious why a high-order process is desired. 

The NEWTON-RAPHSON process, based on 


G(y) = » — Fly)/F'(y) 


is second order, as Hartree points out. In effect, this process draws a tangent 
to the curve Z = F(y) at y = yn, and takes y,4; at the intersection of this 
tangent with the y-axis. The primary source of error is the curvature of F, 
and any operation that reduces the curvature will improve the convergence 
of the iteration. If we write 


H(y) = F(y){F’(y)}-* 


H"(y) = —$F(y){ Fo" (y)(F’)) 7}. 


The function H(y) consequently has the same roots as F, and zero 
curvature at each of them. If we apply the Newton-Raphson process to it 
we get 


we find 


G(y) = y — 2FF'/(2F" — FF"), 


and it is easily verified that this process is third-order. 
As an example, if F(y) = »’ — a we get 


Yati = Yn( Yn? + 3a) /(3yn? + a) 


as a third-order process for computing a*. For a = 10, y = 3, we find 
y1 = 3.16216, ye = 3.162277660168341, which is in error by 4 units in the 
fifteenth figure. 

WILLiAM B. JORDAN 
General Electric Co. 
Schenectady, N. Y. 


on R. HartrEE, “Notes on iterative processes,” Camb. Phil. Soc., Proc., v. 45, 1949, 
p. ; 








QUERIES 


QUERIES 


38. Gauss-SEIDEL Process.—Is there any reason why Gauss’ name 
should be associated with this method of solving linear systems? 
GEORGE E. FORSYTHE 
NBS Institute for Numerical Analysis 
Univ. of California, Los Angeles 


QUERIES—REPLIES 


48. EXPONENTIAL INTEGRALS FOR COMPLEX ARGUMENT (Q 36, v. 5, 
p. 55).—See the tables by E. KrEyszic reviewed in this issue [RMT 909]. 


CORRIGENDA 


, equation (2) for h read k. 
yey (25) for =B* read =aB*. 
1. 22, for (—1%) read (—1)¥. 
, 1. 5, for L read 1. 
, RMT 862, |. 4, for © read — o. 
, RMT 865, |. 4, for x read 2; 1. 6, for x, read x. 





